Universality for Random Tensors 



(N 



o 



- I - 1 

>< 



Razvan Gurau* 
July 11, 2012 



' Abstract 

(N ' 

We prove two universality results for random tensors of arbitrary rank D. We first prove that a 
random tensor whose entries are N D independent, identically distributed, complex random variables 
converges in distribution in the large N limit to the same limit as the distributional limit of a Gaussian 
| tensor model. This generalizes the universality of random matrices to random tensors. 

We then prove a second, stronger, universality result. Under the weaker assumption that the joint 
probability distribution of tensor entries is invariant, assuming that the cumulants of this invariant 
distribution are uniformly bounded, we prove that in the large N limit the tensor again converges in 
Ph distribution to the distributional limit of a Gaussian tensor model. We emphasize that the covariance of 

" the large N Gaussian is not universal, but depends strongly on the details of the joint distribution. 

1 Introduction 



There are two main versions of universality in probability theory. The ordinary version is the central 
limit theorem, stating that the (appropriately rescaled) sum of a large number of independent identically 
distributed (i.i.d.) random variables follows a normal distribution. The second version, or matrix-case, states 
Q\ \ that the statistics of invariant quantities of an N by N random matrix is independent of the details of the 

atomic distribution of the coefficients of the matrix. In the large N limit the random matrix converges in 
distribution to a Gaussian matrix model. In more familiar terms, the eigenvalues density obeys the Wigncr 
semi-circle law under quite general assumptions [TJ [2 [3] . Universality extends to details of the statistics of 
eigenvalues in the large N limit. The spacing of eigenvalues for instance is determined only by the first four 
moments of the distribution of the matrix entries [3] and follows Dyson's sine law J5J[B]. 

In the matrix case the invariant moments are traces of polynomials in the matrix. The limit law can 
be deduced using a Feynman graph representation. In this approach the problem reduces to finding the so- 
called 1 /N expansion for random matrices introduced in [7] . This fixes the correct rescaling of the invariant 
observablcs and their limit distribution. The statistics of the eigenvalue density appears as a clever gauge- 
fixed version of this limit in the particular gauge of diagonal matrices. The apparent non-Gaussian character 
of the Dyson- Wigner law is due to the particular form of the Faddeev-Popov determinant which can be 
computed exactly in this gauge. The resulting Vandermonde determinant governs the eigenvalues repulsion 
hence Dyson's sine law. But universality does not require gauge-fixing. 

Although universality can be established under quite general assumptions, in the matrix case there exist 
invariant probability laws which are not universal [5]. For example any measure which can be written as 
the exponential of the trace of a polynomial in the matrix has a planar but not necessarily Gaussian large 
N limit. A Gaussian matrix can be recovered then via the non commutative central limit theorem. Under 
very general assumptions random matrices become free in the large N limit (this is again a consequence of 
the 1 /N expansion) , and the central limit theorem ensures that the (appropriately rescaled) sum of a large 
number of free matrices converges in distribution to a Gaussian matrix [5J 1101 111) . 

To summarize there arc two ingredients which power both universality and freeness for matrices, namely 
the invariance and the 1 /N expansion. Random matrices encode a theory of random two dimensional surfaces 
and are widely applied in physics for the study of integrable systems, exact critical statistical mechanics, 
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quantum gravity in two dimensions and the list goes on. Matrices generalize in higher dimensions to tensors. 
Introduced in the '90 [12l H3] as tools to study random geometries in dimensions higher than two, random 
tensor models remained an open problem ever since. Although invariant quantities for tensors are well 
known, until recently no 1/N expansion existed for tensors of rank higher than two and no analytic result 
on these models could be established. The lack of results on random tensor is exemplified by the Gaussian 
distribution. One can of course easily write a Gaussian distribution for a random tensor. However its large 
A behavior, that is identifying the appropriate observables (and their scaling) in the large TV limit, has not 
been established prior to this work. 

The situation has drastically changed recently and the necessary ingredients for universality have been 
found for tensors of higher rank, with the discovery of the 1 /A expansion [T?J [T51 [TB] for colored [T71 [H] 
random tensors. The first consequences for statistical mechanics and quantum gravity have been developed, 
see [19j for a general review of this thriving subject. 

In this paper we derive the universality properties associated to this 1/N expansion for a unique complex 
non symmetric tensor. We establish two universality results. The first one is just the straightforward 
generalization of the universality of the Gaussian measure to tensors with entries i.i.d. random variables. 
The second one is more powerful. The natural requirement one should impose on the joint distribution of 
the tensor entries is not independence, but invariance. We show in this paper that if the joint distribution 
of the entries is invariant and its cumulants are uniformly bounded then in the large N limit the random 
tensor converges in distribution to the distributional limit of a Gaussian tensor model. This is in contrast 
with random matrices, and shows in particular that the Gaussian distribution is a more powerful attractor 
for higher rank tensors than it is for matrices. However we emphasize that the covariance of the large A 
Gaussian is not universal and the large A limits of random tensors are rather subtle. The Gaussianity allows 
one only to compute all the large TV correlations in terms of the large A covariance, but the latter has a 
very non trivial dependence on the details of the joint distribution of entries. In particular the perturbed 
Gaussian measures (presented in appendix [A]) lead to a multitude of continuum limits |20j , thus describing 
infinitely refined geometries, dominated by spherical topologies |19| . 

Our results cover tensors of arbitrary rank and lay the foundation for the study of random geometries 
in arbitrary dimensions using random tensors. This study is relevant for critical statistical mechanics, 
intcgrability, quantum gravity and so on in more than two dimensions. 

The proofs of our results rely on a representation of the cumulants of the joint distribution of tensor 
entries by colored graphs. This representation is of course inspired by the Feynman graphs representation 
of perturbed Gaussian measures. However, unlike the former, our representation is completely general and 
applies to all invariant joint distributions of the entries. The precise link between our graphical representation 
and Feynman graphs is detailed in the appendix [AJ Of course, the main challenge is not so much to find 
an appropriate graphical representation, but to compute the contribution of each graph. This requires on 
one hand to find the appropriate scaling of various cumulants with A, and on the other hand a detailed 
combinatorial study of the graphs. If one assumes a uniform scaling of the cumulants (i.e. all cumulants at 
a given order scale with the same power of A, irrespective of the associated graph), the scalings presented in 
this paper arc optimal: tensor distributions which violate them do not admit a large A limit. We comment 
on these scalings and if they can be relaxed in the non uniform case (i.e. when the scaling of a cumulant is 
allowed to depend on the details of the associated graph) in appendix [5] 

One interesting question is to combine our graphical representation with the Connes-Kreimer algebra 
[211 122] of the usual Feynman graphs, as the trace invariant cumulants have the structure of an antipodc 
of a graph Hopf algebra. A second important open question not addressed in this paper is to find a clever 
gauge fixing which would generalize correctly the diagonal condition in the matrix case, and to compute the 
corresponding Faddeev-Popov determinant. This may require to find better "finite- A truncations" of the 
theory (i.e. better cutoffs in the quantum field theory language), and an appropriate generalization of the 
notion of eigenvalues and spectrum for tensors. 

The proofs we present below are combinatorial and rely heavily on the colored graph representation we 
introduce. The plan of the paper is as follows. In section [2] we give the relevant definitions and state our 
two universality theorems. In section [3] we recall the universality for random matrices and its link with the 
1/N expansion. We use this opportunity to introduce at length the colored graph representation for this 
more familiar case. Once familiarized with this representation we present a number of combinatorial results 
concerning colored graphs in the first part of section [U We subsequently use this combinatorial input to 
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prove the two universality results for random tensors in the second part of section 2] Thus the subsections 
I4.1l and r4.2l are mainly review (except lemmas [4] and [5J, while the subsections 14. 31 [4~4l and [4~5l are entirely new 
and contain our main results. In the appendix [Al we give a detailed presentation of the perturbed Gaussian 
measures for random tensors, both in perturbations (subsection IA.1|) for the generic case and at full non 
perturbative level (subsection IA.2[) for a particular example. Both these subsections present new results. 

This paper falls short in many technical points. We do not give a precise definition of infinite tensors, we 
do not propose a generalization of the diagonal gauge of random matrices, we do not detail the sub leading 
corrections in N and so on. All these, and many other, topics need to be thoroughly examined and clarified 
before obtaining a fully fledged theory of random tensors. Our contribution is the derivation of the generic, 
universal behavior of random tensors at leading order, which is the prerequisite for all such studies. 

2 Notations and Main Theorems 

A rank D covariant tensor T n i „d (with n 1 , n 2 , . . . n D 6 {1,...N}) can be seen as a collection of N D 
complex numbers supplemented by the requirement of covariance under base change. We consider tensors 
T with no symmetry property under permutation of their indices transforming under the external tensor 
product of D fundamental representations of U(N). In words, the unitary group acts independently on each 
index of the tensor. The complex conjugate tensor T„i „d is a rank D contravariant tensor 

T' a i„ M D= ^2 U a i n i . . . V a D n DT n i^ n D , Tgi__ _ a o = ^2 U n n n n , . .V n i n if n i_ M n . (1) 
n 1 ...n D n 1 ...n D 

where wc denoted conventionally the indices of the complex conjugated tensor with a bar. We will sometimes 
denote the D-uple of integers n 1 . . . n D by n and assume (unless otherwise specified) D > 3. 

Among the invariants one can build out of T and T we will deal in the sequel exclusively with trace 
invariants. The trace invariants are built by contracting in all possible ways pairs of covariant and con- 
travariant indices in a product of tensor entries. We write such a trace invariant formally as 

Tr(T,T) = ^H<W /„ ...... T n ... , (2) 

where all indices are saturated. By the fundamental theorem of classical invariants of U(N), the trace 
invariants form a basis in the space of invariant polynomials in the tensor entries (see |23j for a direct proof 
relying on averaging over the unitary group) hence in particular the probability distribution of a random 
tensor is encoded in their expectations. 

Remark that a trace invariant has necessarily the same number of T and T and that an index rf is always 
contracted with an index n l . A trace invariant can be represented as a bipartite closed U-colored graph (or 
simply a Z?-colorcd graph). 

Definition 1. A bipartite closed /^-colored graph is a graph B = (V(B),£(B)) with vertex set V(B) 
and edge set £(B) such that: 

• V(B) is bipartite, i.e. there exists a partition of the vertex set V(B) = A(B) U A(B), such that for 
any element I <G £{B), then I = (v,v) with v <G A(B) and v <G A(B). Their cardinalities satisfy 
\V{B)\ = 2\A{B)\ = 2\A{B)\. 

• The edge set is partitioned into D subsets £(B) = U£=i where £ l (B) = {l l = (v, v)} is the subset 
of edges with color i. 

• It is D-regular (all vertices are D-valent) with all edges incident to a given vertex having distinct colors. 

To draw the graph associated to a trace invariant we represent every T n i n n by a white vertex v and 
every T n i n n by a black vertex v. We promote the positions of an index to a color, thus n 1 has color 1, n 2 
has color 2 and so on. The contraction of an index n l on T n \ m „d with an index n % of T fi i fi u is represented 
by a line l % = (v,v) € £' l (B) connecting the vertex v (representing T n i ...„d) with the vertex v (representing 
Tn^...n D )- The lines inherit the color of the index, i, and always connect a black and a white vertex. Some 
examples of trace invariants for rank 3 tensors are represented in figure [1] Every trace invariant can be 
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Figure 1: Graphical representation of trace invariants. 



written as 

D 

Tr B (T,f)=J2€n II T *A' €n = U II 5 KK- (3) 

n,n v,v£V{B) »=1 l i = (v,v)G£ i (B) 

We call the product 8® n encoding the pattern of contraction of the indices the trace invariant operator 
associated to the graph B [24]. The trace invariant associated to a graph B factors over its connected 
components B p . We call a trace invariant whose associated graph is connected a connected trace invariant 
(or a single trace invariant). 

Definition 2. The faces of a D- colored graph B are its connected subgraphs with two colors. We denote 
F^'ii^B) the set of faces with colors i and j ofB, and (B) = \ J : ^> {B) \ their number. The d-bubbles of 
a graph are its connected subgraphs with d colors. 

A colored graph is a cellular complex with cells given by the (i-bubbles. In fact it can be shown that it is 
an abstract simplicial complex, and even more, a simplicial pseudo manifold |17[ 118] , see also appendix |A] 

A random tensor is a collection of N D complex random variables. We consider only even distributions, 
that is the moments of the joint distribution of tensor entries are non zero only if the numbers of T and 
T insertions are equal. We denote the joint moment of 2k tensor entries by niT^jT^. . . .Ta k ,Tf L ^). The 
cumulants of the joint distribution of tensor entries are defined implicitly by 

M 

M(Ttu , Tfii ■ ■ ■ T a k , T n- k ) = II K2fc («) t T "°i '%!-•■]' ( 4 ) 

7r a— 1 

where ir runs over the partitions of the set of 2k points V = {1 . . . k, 1 . . . k} into |a| disjoint bipartite subsets 
Via) = {ai,..., ai, ... } for a — 1, 2, . . . |a|, \a\ < k with cardinality |V(a)| — 2k(a). As the partition 
in a unique set appears only once, the cumulants can be computed in term of the moments. Note that 

We will define a trace invariant distribution as a distribution whose cumulants are trace invariant opera- 
tors. We will allow in this definition trace invariant operators which correspond to disconnected graphs. At 
first sight it might seem rather surprising that according to our definition a cumulant (a connected moment) 
can be expressed as a sum over disconnected graphs. First, the case when the cumulants expand only in 
connected graphs is certainly a particularization of this more general case. Second, and most importantly, 
it is in fact natural to allow disconnected graphs into the expansion of a cumulant in invariants. This is 
clear when dealing with perturbed Gaussian measures in appendix \^ both at the perturbative and at the 
non perturbative level. In perturbations this is seen as follows: moments expand in Feynman graphs, and 
cumulants (connected moments) expand in connected Feynman graphs Q. However the pattern of contrac- 
tion of the tensor indices associated to a Feynman graph Q is encoded in its boundary graph, B = dQ (a 
precise definition of the boundary graph is given section [4.2]) . It turns out that a Feynman graph Q can 
be connected (thus contributing to a cumulant), and have a disconnected boundary graph dQ (as shown 
figure [S]). In order to include the perturbed Gaussian measures one must allow disconnected graphs in the 
expansion of a cumulant. At the non perturbative level this can be seen as a consequence of the invariancc 
of the cumulants under unitary transformations (see the proof of theorem [8]) . The same phenomenon ap- 
pears in the more familiar case of random matrices: at finite TV one obtains contributions to the cumulants 
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corresponding to connected Feynman graphs having two or more external faces ( "multi loop observables" in 
the physics literature). Each external face is a connected component of the boundary graph. However such 
contributions are penalized in the scaling with N. 

We need some more notations. We denote B a generic D-colored graph with 2k(B) vertices labeled 
1, . . . k(B), 1, . . . k(B). We also denote p(B) the number of connected components (labeled B p ) of B, and 
2k(B p ) the number of vertices of the connected component B p . We have Ylp=i M^p) = M^) an d every 
graph B has an associated partition of the vertex set {1, . . . k(B), 1, . . . k(B)} into p{B) disjoint bipartite 
subsets of cardinality 2k(B p ), p = 1, . . . p{B). 

Definition 3. The probability distribution p of the N D complex random variables T,-j is called trace in- 
variant if its cumulants are linear combinations of trace invariant operators, 

p(B) 

K2k [T ni ,f Hl ...T nk ,f H _}= J2 ^^11^' ( 5 ) 

B, k(B)=k p=l 

where the sum runs over all the D- colored graphs B with 2k vertices. 

To compute the joint moments of a trace invariant distribution one has to perform two expansion: first 
the expansion of the joint moments in cumulants and second the expansion of the cumulants themselves in 
graphs. We are interested in the large N behavior of a trace invariant probability measure p^. In order for 
such a limit to exist, the cumulants of pn must scale with N . There are two main cases. Either the scaling 
with N is uniform, that is it insensitive to all but the roughest features of the graph B or it depends on the 
details of B. We will deal in the main body of this paper with the first case, and briefly discuss the second 
case in appendix IB"1 We denote 

m»N) sK{B>N) (6) 



N -2(D-l)k(B)+D-p(B) 



There exists a unique D-colored graph with 2 vertices (all its D lines necessarily connect the two vertices). 
We call it the £>-dipolc and denote it . We call K(B (2 \N) the covariancc of the distribution pjy. 

Definition 4. We say that the trace invariant probability distribution p^ is properly uniformly bounded 

at large N if 

lim K{B {2 \N) = K{B {2) ) < oo , K(B,N) < K{B) , VS ^ B (2) and N large enough . (7) 

N— ¥oo 

We will establish our universality results for properly uniformly bounded distributions. A natural ques- 
tions one can ask at this point is if, in particular examples, proper uniform boundedness is easy to establish. 
This questions is addressed in appendix [A] We first show that uniform boundedness holds perturbatively for 
all perturbed Gaussian measures. Indeed for such measures the cumulants can be expressed as sums over 
Feynman graphs and in appendix I A . 1 1 we show that each graph respects the proper uniform bound. However 
this is not yet a proof: in order to establish proper uniform boundedness of a cumulant one must deal with 
the sum over all Feynman graphs. Sums over graphs are notoriously difficult to control (the perturbative 
series are not summable, but only Borel summable), and promoting a perturbative bound to a bound at the 
full non perturbative level is the object of constructive field theory [25]. We will prove in appendix I A . 2 1 that 
the proper uniform bound on the full resumed cumulants holds for a measure perturbed by a quartic invari- 
ant. The full constructive bounds on cumulants for arbitrary polynomially perturbed Gaussian measures can 
be achieved by an appropriate generalization of the techniques discussed in appendix IA.21 We emphasize 
that once constructive bounds are established they always reproduce the scaling with N of the perturbative 
bounds, hence the perturbative uniform bounds established in appendix I A . 1 1 should hold for the full resumed 
cumulants also in the general case. 

The trace invariance condition of the joint distribution is weaker than the i.i.d. condition. The latter 
can be seen as supplementing the trace invariant operator Ilp^fi ^na by a number of further identifications 
of indices, imposing that all indices of color i in a cumulant are equal (and modifying appropriately the 
scaling with N). These extra identifications decrease the number of independent indices and simplify the 
joint measure. 
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The normalized Gaussian distribution of covariance a 2 for a random tensor is the probability measure 



(T 2 27TI / 

n 

It is characterized by the expectations of the connected (single) trace invariants 

(Tr 6 (T,f)) CT2 =J (n ^" 1<n %f* ) e-^-^^^^^A ^(T.f) . (9) 

n 

It is in fact a non trivial problem to compute the moments of the Gaussian distribution, and we defer it to 
section POl For now we just mention that for any graph B with 2k{B) vertices there exist two non negative 
integers, 0(0) and R(B) such that 

lira N- 1+n< - B U Tr B (T,f)) =a 2k ^R(B). (10) 



N- 

We call fl(B) the convergence order of the invariant B. The normalization in eq. (|8|) is the only normal- 
ization which ensures that the convergence order is positive and, more importantly, for all B, there exists 
an infinite family of invariants (graphs B') such that Cl(B) = fl(B'), see lemma[71 

Definition 5. A random tensor T distributed with the probability measure /in converges in distribution 

to the distributional limit of a Gaussian tensor model of covariance a 2 if the large N limit of the expectation 
of any connected trace invariant equals the large N Gaussian expectation of the invariant 

lim N- 1+n< - B ^ N \Tr B (T,f)} = a 2k ^ R{B) . (11) 

This paper establishes two theorems. The first one simply generalizes the universality of random matrices 
to random tensors: 

Theorem 1 (Universality 1). Let N D i.i.d. random variables Tft, each of covariance a 2 . Then, in the large 
N limit, the tensor = — n^n-Tk converges in distribution to a Gaussian tensor of covariance a 2 . 

The second universality theorem is: 

Theorem 2 (Main Theorem: Universality 2). Let N D random variables whose joint distribution is trace 
invariant and properly uniformly bounded of covariance &(B( 2 \ /Ujv)- Then in the large N limit the tensor 
converges in distribution to a Gaussian tensor of covariance K(B^) = Iinijv-^oo A(B^ 2 \ /Ujv)- 

Universality is thus much stronger for random tensors than it is for random matrices. For the latter 
universality can be established if, for instance, the distribution /ijy is i.i.d, but one achieves various non 
Gaussian large N limits [8] for trace invariant measures. The limit eigenvalue distributions can be evaluated 
and it is different from the usual semicircle law (multi cut solutions and so on). A set of matrices whose joint 
distribution is trace invariant become free in the large N limit. Random tensors exhibit a more powerful 
universality property: properly uniformly bounded trace invariant distributions become Gaussian in the 
large N limit. However note that the large N covariance K(B^ 2 ') strongly depends on the details of the joint 
distribution at finite N. For the case of perturbed Gaussian measures the large N covariance is a sum over 
an infinite family of Feynman graphs and exhibits various multicritical behaviors [20j . 

Before proceeding we fix some notations. From now on B will always designate the invariant whose 
expectation we evaluate. As we deal only with connected (single trace) invariants, B will always be a 
connected D colored graph. The graphs B(a) arise from the expansion of cumulants into trace invariant 
operators. They are not connected. Their connected components are labeled B p (a). 

When evaluating expectations of observables we will introduce D + 1 colored graphs (definition [T] with 
D replaced by D + 1). We will call the new color 0. We will use Q as a dustbin notation for connected D + 1 
colored graphs. The lines of the new color 0, denoted 1° E £ (Q), play a special role and will be represented 
as dashed lines. 
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3 Random Matrices 



We we will first detail the case of random matrices. This serves both as motivation and as an opportunity 
to introduce the appropriate tools for the study of random tensors. 

All connected bi-colored graphs with 2k vertices are cycles with alternating colors (which we denote B) 
The associated trace invariants write 

2 

S B -=Y\ 17 6 t-i 

w nn n v n v 

i=i i*=(«,«)e5*(B) 

Tr 6 (A, A) = £ C II = Tr[(AAt) fc ] , (12) 

n,n v,veV(B) 

Any invariant function of a generic (i.e. not necessarily hcrmitian) matrix can be evaluated starting from 
these trace invariants, as they fix the spectral measure of AA^ . 



Gaussian distribution of a random matrix. The Gaussian distribution of a non hermitian random N x JV 
matrix A of covariance 1 is the probability measure 



J2 A„i„2 <5„i fl i j„ 2ft2 A 5 i 



J[ (jV ^" 1 ^" 1 " 2 ) , (13) 

(n\n 2 ) 



where the product is taken over all the (complex) entries A„i„2. Note that the exponent can alternatively be 
written in the more familiar form iVTr(AA^). The Gaussian distribution is characterized by its expectations 
in the large TV limit, 

lim AT- 1 /Tr[(AtA) fe ]\ = — !— f 2 , ) , (14) 



JV->oo \ 1 1 1 k + l\k 

It is instructive to prove this. We represent the trace invariant as a colored cycle B with 2k vertices 

Tr[(AtA) fe ])=^C( J] Aa.Ag,). (15) 

n,n v,v£V(B) 

The Gaussian expectation of a product of matrix entries is a sum over pairings (Wick contractions in the 
physics language) of products of covariances. If two matrix entries are paired by a covariance we connect 
them by a dashed line (to which we associate by convention the color 0). A pairing is then represented as a 
(Feynman) graph Q. 

Definition 6. We call a graph with 3 colors Q a covering graph of B if Q reduces to B by deleting the 
lines of color 0, G\£ {G) = B. 

The contraction of two entries A„i„2 and A R i fi 2 with the Gaussian measure (| 13[) comes to replacing them 
by jf5 n ini5 n 2 n 2, hence each line of color 0, 1° — (v, v) G £°(G), will bring a factor jj5 n i n i S n 2 n 2, 

The graph of the invariant B has two colors 1 and 2, while a covering graph Q has three colors: 1, 2 and 
the extra color of the dashed lines. An example of a covering graph Q contributing to the expectation of 
Tr^A'fA) 3 ] is presented in figure [5] 

The expectation of B becomes a sum over all covering graphs Q 

Tr[(AtA) fc ])=£(n II S K*i) E II ^ 6 <^Mn%, (16) 

n,n i=ll*=(v,v)EE*(B) G,G\£°(G)=B l = (v,v)<££°(Q) 

and, as the lines of color 1 and 2 of any such Q are in fact the lines of color 1 and 2 of B 

(Tr[(AtA) fc ])= E E (IT II <^n|)( II ^nlnlKmi)- (17) 

G,G\£°{G)=B n,n t=l p^v^&S^g) l° = {v,v)e£°(.G) 
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Figure 2: A covering graph Q of an observable B. 



To evaluate the contribution of a graph Q one must evaluate the number of independent sums over the 
matrix indices n, n. The Kronecker S compose along the faces (bi-colored circuits) of colors 01 and 02 and 
yield an independent free sum for each such face. As we have exactly k lines of color we get 

(Tr[(AtA) fe ])= J2 ^ F ° 1(e) + f ° 2(e) ■ (18) 
G,G\£°(G)=B 

Note that the face 12 corresponding to the circuit B with colors 12 (hence to the observable itself) does not 
bring any sum. The graph Q has 2k vertices (k black and k white), 3k lines (k dashed lines of color and k 
solid lines for each of the colors 1 and 2) and faces (F 01 (5) + F 02 (Q) representing free sums and F 12 (Q) = 1 
with no sum) . The Euler character of Q is 

2k -3k + F Q1 (G) + F O2 (0) + 1 = 2 - 2g(G) => -1 - k + F al (g) + F 02 (g) = -2g(g) . (19) 

It follows that in the large A limit only graphs g of genus g(g) =0 contribute. We call such graphs minimal 
covering graphs of B. Equivalently they can be seen as the covering graphs of B with maximal number of 
faces F 01 (g)+F 02 (g). Thus 

lim Ar 1 /Tr(A t A) fc \ = R k , (20) 

where Rk counts the number of minimal (planar) covering graphs g, g \ £°(g) = B. It is easy to see that 
Ri = 1 and Rk+i = 2p=o RpR-k-p, thus Rf. = ^-j- ( 2 / , fe ) , i.e. R k are the Catalan numbers. The normalization 
of the Gaussian is canonical, and not a matter of choice: any other normalization leads either to infinite or 
to zero expectations in the large limit. 



3.1 Universality for Random Matrices 

In order to introduce the ideas we will use later to prove the universality properties of random higher rank 
tensors we present below the classical universality of random matrices using this graphical representation. 

Theorem 3. Let M be a matrix with entries i.i.d. complex random variables with centered distributions of 
unit covariance. In the large N limit, the matrix M = converges in distribution to a random matrix 

distributed on a Gaussian. 

Proof: Non hcrmitian matrices whose entries are i.i.d. complex random variables are called random covari- 
ance matrices [8]. The moments of the matrix M write 



lim 

N— ¥00 



lim 



AT ^ [ Tr 
1 v , 



= lim — — r 

N^oo N 1 + k 



UN 



Tr(MM*) 1 



nil 



[ n 



M*M S _ 



tVDSV(E) 



(21) 
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where we denoted k t the product of cumulants associated to the partition ir. As the entries are independent, 
the only non zero cumulants are K 2q (Mij) q (Mij) q . Like in the Gaussian case, each cumulant will introduce 

constraints on the number of independent sums. We slightly extend our graphical representation. If two 
matrix entries are connected by a two point cumulant we connect them, as in the Gaussian case, by a 
dashed line of color 0. If four (or more) matrix entries are connected by a cumulant, all the four (or more) 
matrix elements have the same indices. We will employ a simple trick to represent such cumulants, namely 
we will connect the matrix entries two by two (a M and a M) by dashed lines of color and keep in 
mind that the indices are further identified. The pairing is not canonical, and in order to control the sub 
leading contributions one needs to improve this graphical representation and track carefully the higher order 
cumulants. However at leading order we just need a rough estimate of the number of independent sums in 
an observable and a non canonical pairing suffices. 

The graphs Q we obtain are covering graphs of B, Q \£°(Q) = B. We have (at most) an independent 
sum over an index corresponding to the faces 01 and 02 (potentially less if several dashed lines correspond 
to a higher order cumulant). In the large N limit only planar graphs (minimal covering graphs) contribute. 
Furthermore, if such a planar graph corresponds to a factorization with a fourth (or higher) order cumulants, 
some of the faces 01 and 02 arc further identified (as a pair of distinct lines of color on a planar graph with 
a unique face 12 can never share both faces 01 and 02), hence the number of independent sums is strictly 
smaller than F 01 (Q) + F 02 (Q) in this case. It follows that the only surviving contributions in the large N 
correspond to planar graphs in which all dashed lines come from a second order cumulant 



lim — u 

JV-s-oo N 



Tr 



lim 

JV->-cx 



g,Q\£°(g)=B 



K 2 



(22) 



where we used the fact that the covariance of the atomic distribution is one. 



□ 



In the case of matrices we have another clever set of observables, the eigenvalues of the matrix 
which are non-polynomial functions of the generators. Passing to this set of variables is analog to writing the 
theory in a particular gauge and the corresponding Faddev-Popov determinant results from the integration 
over the unitary group with the Haar measure. The result is the well known Vandcrmonde polynomial. 

We now relax the requirement of independence and require only trace invariance of the joint distribution 
of the entries. Thus in eq. ([2"TT) 
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Ar->oo TV 



Tr 
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x) N 
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(23) 



one substitutes for each set in the partition 7r the properly uniformly bounded trace invariant cumulants of 
cq. © and © 



K 2kia) [M ni ,M Kl ...M n 



p(b(«)) 

£ N -*W°))+*-p( B < a VK(B(a),N) JJ S& a) . (24) 



B(a), k{B(a))=k{a) 



p=l 



The index a = 1, . . . \a\ tracks the cumulant K 2 fc( a ) appearing in the expansion of the joint moment. The 
index p = 1, . . .p(B(ay) labels (at fixed B(a)) the connected components B p (a) in the expansion of K 2 k(a) 
in trace invariants. 

When evaluating the expectation of a trace observables, the sum over partitions ir becomes a sum over 
graphs Q. The graph Q representing a term in the sum is constructed as follows. First one draws the 
observable B and an invariant B(a) (with connected components B p {a)) for each K 2 k( a ) f° r a = 1, • • • 



Note that Yfp=i k(jB p (a)) = k(a) and X)lT=i k(a) — k. As a matter of convention we flip all the black 
and white vertices of B. Note that in this graphical representation all the original vertices of B are doubled: 
every vertex appears once in B and once in some B p (a). We connect every vertex representing a matrix entry 
M in B with the vertex representing the same matrix entry M in the corresponding B p {a) by a fictitious 
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dashed line of color 0. We explicit the label a of the TD-colored graphs B p (a). Some example are presented 
in figure [3] 

We thus construct a closed connected graph Q having three colors, 0, 1 and 2. As we nipped the black 
and white vertices on B, all lines of color in Q will connect a black and a white vertex. We call a graph 
built in this way a doubled graph. The sums over partitions 7r and invariants B(a) in equations (|2"5|) 
and (|24[) becomes a sum over all doubled graphs Q one can build starting from B which we denote Q Z> B. 
Starting from a given Q one readily identifies B,B p (a) and p(B(a)): the observable B is the subgraph with 
colors 1, . . . D of Q having no label a, all the other subgraphs with colors 1, . . . D of Q represent the various 

£? p (a)'s, that is Q \ £°(Q) = B U^=i (y P p=i Bp( a )J s an( i p(B(a)) is the number of connected components 
of Q sharing the same label a. 

This graphical representation applies to all trace invariant measures. We will see in appendix [5] the 
precise relation between the usual Feynman graphs for perturbed Gaussian measures and these doubled 
graphs, but we warn the reader that this relation is more subtle than it might appear at first sight. 

Some doubled graphs contributing to the observable Tr[(MM^) 3 ] are given in figure^ The face 12 
associated to B is the one with six vertices, while the faces 12 with four and two vertices correspond to 
various B p (a). We have also identified on the drawings the various cumulants a to which each connected 
component B p (a) belongs. Thus on the left hand side of figure [3] we represented a contribution form two 
cumulants. The first one is a two point cumulant k(B(l)j = 1, and the second one is a four point cumulant 
fc(£>(2)) = 2. The invariant for the first cumulant has a connected component = 1 with two vertices 

k(Bi(l)) = 1. The invariant for the second cumulant has also one connected component p(S(2)) = 1 but this 
time with four vertices fc(Z3i(2)) = 2. On the right of figure[3]we presented a contribution coming from the 
same two cumulants, A;(i3(l)) = 1, fc(£>(2)) = 2. The invariant for the first cumulant has again a connected 
component p(B(l)) = 1 with two vertices k(B\{V)) = 1. But this time the invariant for the second cumulant 
has two connected components p(2?(2)) = 2, each with two vertices k(Bi{2)} = 1, fc(S 2 (2)) = 1. 




Figure 3: Doubled graphs contributing to an observable. 



To evaluate the contribution of a graph Q to the expectation of an observable one must remember that 
we first divide the 2k points among |a| cumulants, and subsequently the 2k(a) points in every cumulant are 
subdivided into p(B(a)) connected graphs B p (a). As the lines of color connect two copies of the same 
vertex, the indices of their end points are identical, hence each I = (v,v) £ £°(G) contributes 5 n i n i5 n 2 n 2. 

The expectation of an invariant observable becomes 



J_ >p jv EL°il(-2fc(e(a)) + 2-p(s( a ))^ ' 

GDB, e\£°(e)=BuL°i 1 (u^"') B„(a)) 



gdb, g\e°(g)= 

\a\ p(B(a)) 



X 

n.n 



EK*n n s nn) n <w^- ^) 

a=l p=l l° = (v,v)££°{g) 

The total operator (s^ n Yl^=i Tl^i ^n^"' > ) explains our representation in doubled graphs: one must 
keep track of the observable B, the cumulants K2fc(a) and the graphs B p (a) in order to compute the contribu- 



te 



tion of a term to the expectation of the observable. In particular this requires the doubling of the vertices. 
Substituting the trace invariant operators 8f lfl and eq. (|25p becomes 

I £ ^(-^N^N) JV) 

CDS, e\£0(e)=Bu^i 1 (u^ (Q) ' , B p (a)) 
2 

Z- = (^0)6£^BuL Q i 1 (u p i 1 ; B P ( Q ))J 

and noting again that the lines of colors 1 and 2 of B U^Li {^p^i ^p( a )) arc exactly the lines of color 1 
and 2 in Q, we see again that the Kronecker S's compose along the faces of colors 01 and 02 of G, thus 

l^(Tr(MtM) fc ) = J2 Ar- 1 - 2fc+ ^^ p ( B («))+ i?01 ^+^ 02 (S)-2EL a i 1 (p(B(a ! ))-i) f[ K (B{a),N). (27) 

G,GDB a=l 

The doubled graph Q has 4fc vertices, 2k coming from B and 2k coming from all the B p (a). It has 

1 + X)lr=ii p{&( a )) f aces 12, one associated to the observable B, and one for each B p (a). Furthermore it has 
2k lines of color 0, 2k lines of color 1 and 2k lines of color 2. The Euler character of Q is 

M 

4fc - 6k + 1 + J2 P( S («)) + f01 (S) + F 02 {g) = 2 - 2g(Q) , (28) 

Q = l 

hence the global scaling with N of a term is iV _29(e)_2 ^" 1 ( p ( e(Q) )~ 1 ) . It follows that Q contributes to 
expectation of an observable in the large N limit if it is planar and each cumulant K 2 k{ a ) contributes exactly 
one connected invariant p(B(a)) = 1. The second condition is easy to understand for perturbed Gaussian 
measures. As previously stated the disconnected invariants B(a) correspond to Fcynman graphs having more 
than one external face. Reconnecting the external lines on such a cumulant on the observable B leads to 
non planar Feynman graphs, in spite of the fact that the associated doubled graph (which only sees the 
boundary of the Feynman graph contributing to the cumulant) is planar. This emphasizes the non trivial 
relation between Feynman graphs and doubled graphs. 

The planar graphs contributing to the large N limit possess cumulants of orders between 2 and 2k (each 
cumulant contributing only when its associated invariant is connected), hence the large N distribution of 
M is not Gaussian. The restriction of trace invariant measures for matrices to planar graphs has a different 
effect: one can easily show that matrices distributed according to such measures become free in the large N 
limit. This is particularly transparent in the combinatorial formulation of free probability theory of [261 127j . 
In the large iV limit only the free cumulants (defined by restricting the sum in eq. (|4]) to non crossing 
partitions) survive, and one can show that (in the large TV limit) the mixed free cumulants of a collection of 
matrices cancel. As one only deals with the N — > oo limit, the free cumulants are automatically associated 
to connected invariants. One example of a random matrix model whose measure is not trace invariant is the 
Grosse Wulkenhaar model |55] which is only almost trace invariant. 



4 Random Tensors 

We now go to the core of our paper and the proofs of the two theorems. We start by an account of properties 
of D and D + 1 colored graphs we will use in the sequel. Most of the lemmas we present in subsections 14.11 
and 14.21 can be found in [T^l HH1 HH] • The rest of this section is new. 
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4.1 D + 1-colored Graphs 

The connected (single trace) observables of tensor models are represented by connected .D-colored graphs B. 
Their expectations are evaluated in terms of D + 1-colored graphs Q, having an extra color 0. We will use 
the shorthand notation = {1, . . . D}. 

Consider a connected D + 1 colored graph Q . To simplify notations we will drop in this subsection as 
much as possible Q from our notations. Thus the sets of vertices, edges and faces of colors ij (definition [2]) of 
Q are denoted V, £ and J 7 ^ 1 ^' . Furthermore we denote T = U^jJ 7 ^'^ and F = \1F\. We define the jackets 
[13 HH1 H9] of the D + 1-colored graph Q. 

Definition 7. A colored jacket J is a 2-subcomplex of Q, labeled by a (D + 1)- cycle r, such that: 

• J and Q have identical vertex sets, V(S) = V; 

• J and Q have identical edge sets, £(J) = £ : 

• the face set of J is a subset of the face set of Q: T(J) = U^L ^°^ T . 

For example the jacket associated to the cycle (0,1,2... D) contains the faces (01)(12)(23) . . . (DO). It is 
evident that J and Q have the same connectivity. A given jacket is independent of the overall orientation 
of the cycle, meaning that the jackets are in one-to-two correspondence with (D + l)-cyclcs. Therefore, the 
number of independent jackets is Dl/2 and the number of jackets containing a given face is (D — 1)!0 

The jacket has the structure of a ribbon graph, |30j . as each edge of J lies on the boundary of two of its 
faces. A ribbon line that separates the two faces, (t~ 1 (i), i) and (i,r(i)) inherits the color i of the line in Q. 
Ribbon graphs are well-known to correspond to Riemann surfaces [30] . and so the same holds for jackets. 
Given this, we can compute the Euler character of the jacket, x(J) = \F(J)\ — \£(J)\ + \V(J~)\ = 2 — 2g(J), 
where g(J) is the genus of the jacket 

Definition 8. The convergence degree (or simply degree,) of a graph Q is lu(Q) = ~}2jg(J), where the 
sum runs over all the D\/2 distinct jackets J of Q. The degree is a nonnegative integer. 

Consider a jacket J of a (D + 1) colored graph Q with 2k = \V\ vertices. The number of vertices and 
lines of J are: IV^)! = |V| = 2k and |£(J0| = |£| = (D + l)k, respectively. Hence, the number of faces of 
J is 1^(^7)1 = (D — l)k + 2 — 2g(J). Taking into account that Q has ^D\ jackets and each face belongs to 
(D — 1)1 jackets we obtain 

F = I-FI = ^ E |H7)| = ^p±k + D - j^ U (S) ■ (29) 



This equation is crucial in establishing the universality results in the large N limit of random tensor models. 
Of course the same equation holds (replacing D by D — 1) for Z?-colored connected graphs. Note that, as F 
is an integer, ui(G) is a multiple of ny^jj\ ■ 

We now consider the D-bubbles of Q with colors (i.e. the connected subgraphs of Q with lines of colors 
1, 2 ... D). We denote them B^) . As they are .D-colored graphs, they also possess jackets, which we denote 
by J^y It is rather elementary to construct the jackets of the bubbles from the jackets of the graph J 
[HI HH [19] . Let us construct the ribbon graph J° consisting of vertex, edge and face sets: 

V(J d ) = V(J) = V, £(J d ) = £(J) \ £°(J) =S\£°, 

that is having all the vertices of Q, all the lines of Q of colors different from and some faces of Q. For 
instance, for the jacket corresponding to (0, 1, . . . D) the ribbon graph J° has faces (12) ...(£)— ID) and 



1 It is however sometimes more transparent to over count the distinct jackets by a factor of two associating them one to one 
with cycles. For example, on can count that from the D! cycles of D + 1 colors, (D — 1)! will contain the pair ij and (D — 1)! 
the pair ji. 

2 A moment of reflection reveals that the jackets necessarily represent orientable surfaces. 
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(Dl). Given that the face set of J is specified by a (D + l)-cycle r, the first thing to notice is that the 
face set of J° is specified by a I?-cycle obtained from r by deleting the color 0. The ribbon graph J° is the 
union of several connected components, J^y Each J®^ is a jacket of a £>-bubble £>( M )- Conversely, every 
jacket of 23( M ) is obtained from exactly D jackets of Cid. 

Lemma 1. Let Q be a closed connected D + 1 colored graph and B(u\ its D-bubbles with colors 0. Then 

w(£)>D $>(!?<,.))• (31) 

As are in one-to-one correspondence with disjoint subgraphs of J we have gj > g^-g . As every 

jacket is obtained as subgraph of exactly D distinct jackets J , summing over all the jackets of Q proves 
the lemma (see [H] for more details). 

Of particular importance in the sequel are the graphs Q of degree zero, ui(G) = 0. They have been 
extensively discussed in [55]. In D > 3, the D + 1 colored graphs with degree zero have a very simple 
structure. A counting argument proves that such a graph must have at least a face with exactly two vertices. 
As all the jackets must be planar this in turn implies that the graph contains two vertices separated by 
exactly D lines. Albeit simple, the proof of the second statement is somewhat convoluted. 

For 2+1 colored graphs the degree equals the genus of the graph, hence the graphs of degree are the 
planar graphs. For D > 3, the D + 1 colored graphs of degree zero are called melonic. 

Lemma 2. Let D > 3. If Q is a closed connected D + 1 colored graph of degree zero then Q has a face with 
exactly two vertices. 

Proof: Since Q is of degree zero it has F = DI ' D 2 1 - > k + D faces, from equation (|29l) . Denote F s the number 
of faces with 2s vertices (every face must have an even number of vertices). Then 

F 1+ F 2+ Y J F s = D{D 2 1) k + D. (32) 

s>3 

Let 2fc^ be the number of vertices of the /Lt'th face with colors if. We count the total number of vertices by 
summing the numbers of vertices per face iK - k 1 ?-. = Fi + 2F 2 + J2 S >3 s ^* = DI " D ) +1 ' > k (as each vertex 
contributes to D(D + l)/2 faces). Substituting F 2 from (f3"2"| we get 

F 1 = 2D + J^is - 2)Fs + D{D 2 3) k . (33) 

s>3 

Notice that on the right hand side, the first two terms yield a strictly positive contribution for any D > 2, 
whereas the third term changes sign when D = 3. 

□ 

This lemma explicitly breaks when D = 2: there exist planar graphs having no face with exactly two 
vertices. This is the deep origin of the fact that trace invariant measures can lead to non Gaussian matrices, 
but (as we will prove below) necessarily lead to Gaussian tensors in the large N limit. 

Lemma 3. If D > 3 and Q is a closed connected D + 1 colored graph of degree zero, then it contains a 
D-bubble (i.e. subgraph with D colors) with exactly two vertices. 

We emphasize that the D lines of the D-bubble with two vertices can have any colors, 1, . . . ,D but also 
0,2,...L> or 0, 1,3, ...D, etc. 

Proof: From the previous lemma Q has a face (say of colors if) with exactly two vertices (say v and v). If, 
for all q, a unique line of color q connects v and v we conclude. If the two lines of color q are different, v £ l\, 

3 A jacket iT^j of ^(m) ' s s P ec ifi e d by a D-cycle (missing the color 0). On can insert the color anywhere along the cycle 
and thus get D independent (D + l)-cycles. 
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v <G l 2 we consider the jacket J = (. . . iqj . . . ). It contains the faces (iq) and (gj). As Q is of degree zero, J is 
planar. As Zf and l\ separate the same two faces iq and qj, the graph J' obtained from J by deleting them 
has two lines less, but the same number of faces as J . The Euler character of J' is x(i/') = x(i7) + 2 = 4, 
hence J' has two planar connected components. Then l\ and l\ separate a two point graph having at 
least two vertices less than Q. Note that is not a closed D + 1 colored graph, as it has two vertices 
which are not touched by lines of color q (the vertices on which l\ and l 2 were hooked). 

We now take aside . It can be transformed into a genuine closed D + 1 colored graph by adding a 
line l\ 2 obtained by reconnecting If and l 2 . Consider any jacket Jg^y of . It is a planar ribbon graph 
with one external face. Adding l\ 2 leads to a ribbon graph J7g(g) having one more line and one more face 
(the external face of Jg( q y is divided into two faces by the new line l\ 2 ). It follows that Jg(q) is planar. The 
ribbon graph JgU) is one of the jackets of Q^ q \ hence oj(Q^) = 0. 

Note that one can not naively iterate the argument, as the graph has a line, l\ 2 , which does not 
belong to Q. However, Q( q > has a face of colors i'j' with exactly two vertices v',v'. Again for all q' wc 
consider the graph Q^ q ' q > obtained from Q^> by reconnecting the two lines of color q' containing v' , v' into a 
unique line l\ 2 . The line l\ 2 belongs to only one of these graphs, for a fixed q' . We then chose another one, 
say Q( q ' q ) to iterate (if for all q" ^ q' the two vertices are connected by a unique line we obtained a bubble 
of Q with exactly two vertices and conclude). As l\ 2 is not a line of Q^ q ' but l\ 2 is a line of g (q ' q "\ all but 
one of the lines of Q^ q,q - 1 belong to Q. We iterate until we reach a graph Q^ q,q '•" ) with exactly two vertices 
connected by D + 1 lines. Out of them D are lines of Q and form a D bubble. 

□ 



4.1.1 Melons 

We call two vertices separated by D lines in a graph with D + 1 colors a melon (or an internal -D-dipole, not 
to be confused with the D-dipolc B 1 - 2 * 1 ). Wc emphasize that a melon can have external legs of any color 0, 1 
up to D. The D internal lines of a melon with external lines of color i have colors 0,1, . . .i — 1,1 + i, . . . D. 
Replacing a melon by a line corresponding to its external legs wc obtain a graph of degree zercQ having 
two vertices less (and D ^ D ^ 1 ^ less faces). Iterating, one reduces a graph of degree zero to a graph with 
exactly two vertices connected by D + 1 lines. Conversely all graphs of degree zero can be built by arbitrary 
insertions of melons on lines. The graphs of degree zero are then in one to one correspondence to colored 
rooted D + 1-ary trees [24l [29] . 

First order. The lowest order graph consists in two vertices connected by D + 1 lines. We represent 
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Figure 4: The first order melonic graph and its corresponding rooted tree. 



this graph by the tree with one vertex decorated with D + 1 leaves. A leaf is a vertex with only one incident 
edge. The D + 1 leaves correspond to all the edges incident to the black vertex v. (of course they are all also 
incident to the white vertex v). The leaves inherit the colors. This first vertex is called the root vertex (and 
is marked R) . We consider all edges incident at the black vertex to be active. The leafs of the tree inherit 
this activity. See Figure 0] for an illustration. 

Second order. At second order, D + 1 graphs contribute. They arise from inserting a melon (that is 
two vertices connected by D lines) on any of the D + 1 active lines of the first order graph. Say, we insert 
the new melon on the active edge of color 1. With respect to the new melon, all edges incident at its black 
vertex are deemed active, while the exterior edge (of color 1) incident at its white vertex is deemed inactive. 
This graph corresponds to a tree obtained from the first order tree by connecting its leaf of color 1 to a new 
(D + 2)-valcnt vertex. This new vertex has D + 1 leaves, one of each color. The root and the new vertex 
are joined by tree line of color 1. The leaves correspond to the active lines (either of the root or on the 

4 Every jacket has two less vertices, D + 1 less lines and D — 1 less faces, hence its genus does not change. 
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Figure 5: A second order melonic graph and its corresponding tree. 



new melon). We presented this in figure [5] The inactive line of the graph (represented in bold in figure [5]) 
corresponds to the tree line. All the active lines of the graph correspond to the leaves of the tree. 

Order k + 1. We obtain the graphs at order k + 1 by inserting a melon on any of the active lines of 
a graph at order k. Once again, with respect to the new melon, all edges incident to its black vertex are 
deemed active, while the exterior edge incident to its white vertex is deemed inactive. In terms of the trees, 
we represent this insertion by connecting a (D + 2)-valcnt vertex, with D + 1 active leaves, to one of the 
active leaves of a tree at order k. The new tree line inherits the color of this leaf. 

The 2k vertices of the graph are in two to one correspondence to the k vertices of the tree. The (D + l)k 
lines of the graph are in one to one correspondence to the (k — 1) lines and Dk + 1 leaves of the tree. The 
tree associated to a graph is a colored version of a Gallavotti-Nicolo tree [3T] . 

If a graph is a (D + l)-colored melonic graph, all its subgraphs with D-colors (D-bubbles) are melonic. 
This is easy to see from the construction algorithm. Moreover, the D-aiy trees of the D-bubbles with colors 
0, are trivially obtained from the (D + l)-ary tree of the graph Q by deleting all lines and leaves of color 
0. 

We will use in the sequel the following two lemmas. 

Lemma 4. Let a melonic D-colored graph B. Then there exists a unique melonic D + 1 colored graph Q 
with the same number of vertices which reduces to B by deleting all the lines of color 0. 

The unique D + 1-ary tree Tg with k vertices which reduces to a given D-ary tree 7s with k vertices by 
deleting all the tree lines and leaves of color is the tree 7b decorated by a leaf of color on each of its 
vertices. 

Lemma 5. Let a melonic D-colored graph B with 2k vertices. Then there exists a unique melonic D + 1 
colored graph Q with 4fc vertices which reduces to B by deleting all the lines color 0, such that no two vertices 
of B are connected (when seen as vertices in Q) by a line of color 0. 

As no two vertices of B are connected (in Q) by a line of color zero, if follows that none of the tree vertices 
of the tree 7b associated to B (when seen as a subtree of the tree Tg associated to Q) has a leaf of color 0. 
Therefore all the vertices in 7b must be connected in Tg to another vertex by a tree line of color 0. The tree 
Tg, obtained from 7b by decorating each vertex with a line of color (and a new end vertex), is unique and 
so is its associated graph Q with 4fc vertices. 

4.2 Open graphs and the boundary graph 

We have discussed so far closed connected D + 1 colored graphs. We will now present open D + 1 colored 
graphs, that is graphs having some external lines. 

Definition 9. A bipartite open D + 1-colored graph is a graph Q = (y{Q),E{Q)^j with vertex set V(Q) and 
edge set £(G) such that: 

• V{G) is bipartite, i.e. there exists a partition of the vertex set V{Q) = A(G) U A{G), such that for 
any element I € £{G). then I = (v,v) with v G A(G) and v G -A(G). Their cardinalities satisfy 
\V(G)\ = 2\A(G)\ = 2\A(G)\. 
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• The positive (negative) vertices are of two types, internal vertices and external vertices, A(G) = 
•Aint(G) U^4e X t(0) } A(G) — Ai n t(G) UAcxtiG)- The internal vertices areD + l valent while the external 
vertices are 1-valent. 

• The edge set is partitioned into D subsets £{B) = (J^ £ l (B), where £ l (B) = {l l = (v,v)} is the subset 
of edges with color i. Furthermore the set of edges of color is partitioned into internal and external 
lines of color 0, £°(G) = £i n t(S) £o X t(G) , such that the internal lines connect two internal vertices and 
the external lines connect an external and an internal verte^. All the lines of color i are internal. 

• The lines incident to a D + 1 valent internal vertex have distinct colors, while the line incident to an 
external 1-valent vertex has color 0. 

Some examples of open 3 + 1 colored graphs are presented on the left in figure [6] Both graphs have four 
external lines and four external vertices. 
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Figure 6: Open graphs and their boundary graphs. 



The faces with colors Oi of an open D + 1 colored graph (defined still as subgraphs with lines of colors 
and i) fall in two categories: either they are internal faces, denoted jf^t (l^nt^l = ^int) i- c - they contain 
only internal lines, or they are external faces, denoted ^Vxt* (l^ext^l = -^ext) i- e - they contain external 
lines of color 0. 

The external faces / € ^ext necessarily start and end on two external vertices u and u, f = (u,u). For 
every graph G we build the boundary graph dG having a vertex u (resp. u) for every external vertex of 
G and a line of color i joining a u and a u for every external face / = (it, u) G ^ext °f On ^ ne right in 
figure [5] we represented the boundary graphs dG of the two graphs G- The boundary graph is a D colored 
graph and represents a tensor invariant, thus TT f _, r (o,«) <5„« n i_ = 5^ fl . 

Note that, as it is that case in the second example, in spite of the fact that G itself is connected, the 
boundary graph dG can be disconnected. 

4.3 Gaussian Distribution for Tensors 

We now compute the large N trace invariant moments of the Gaussian distribution for a random tensor 

(Tr B (T,f)) = y g g^ )e-^^ Tr e (T,f) , (34) 

ft 

with the connected trace invariant operators 

Tr B (T,f) = ^^ fi J] T Hv % v , 5l- L = TJ 5 <R% , (35) 

indexed by connected graphs B with colors 1 . . .D having 2k(B) vertices (and Dk(B) lines). Assume a = 1. 
The number of faces of the D-colored graph associated to the observablcs computes from cq. |29|) in terms 

5 Or two external vertices. 
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of its degree 

£ F ij {B ) = {D l) l D 2) k(B) + (D-1)- j^L- u{B ) . (36) 



l<i<j 

The Gaussian expectation is a sum over contractions. As in the matrix case, we represent two tensors 
connected by a covariance as a dashed line to which we assign the color 0. We denote the full graph, including 
the color by Q . An observable is a sum over graphs Q which restrict to B by erasing the dashed lines of 
color 0. We already encounter such graphs in the case of matrices. 

Definition 10. A D + 1 colored graph Q is called a covering graph of B if it reduces to B by erasing the 
lines of color 0,Q\ E°{0) = B. 

Every face of colors Oi in Q brings a free sum, hence a factor N. Every dashed line generated by the 
covariance brings a factor N n-i ■ The moments of the Gaussian write 

Tr B (T,T)\ = J2 N~ k ^ D ~»N^ F ° z ^ 

g, g\s°{g)=B 

N~ HB ^ D ~ 1)+ ^°<'<i Fii( - 9) ~^<'<i Fi:i( - s) , (37) 

9, S\£°(S)=B 

Note that X)o<i<j F^(Q) is the total number of faces of the graph Q, while J2i<i<j F 1 '' {Q) = Ei<i<j F % i(B) 
is the number of faces of B. Also we have k(Q) = k(B). Using eq. ([29]) and (|36|) . we get 

Tr B (T,T)\ = N^T^W^+T^W^ . (38) 

g, g\£°(g)=b 

As both (]j1iy ^{G) and ^ D 2 _ 2 y U>(B) are integers, the scaling with A of a graph Q contributing to the 
expectation of a trace invariant is always an integer. By Lcmma[TJ lu{Q) > Du{B), thus 



2 2 



^Q) + 7T^ o\T w (^) = 1 " ^(S) - TTTTt o\T "(S) - Dlo(B) < 1 - — U (G) . (39) 



(D-l)\ y ' (D-2)\ y ' D\ w/ D{D-2)\ 



4.3.1 Melonic observables 



2 



D\ 



From cq. (|38|) and (p9|) it follows that in the large N limit the expectation of an observable scales at most 
like N, and it scales like N only if there exists a melonic graph Q which restricts to B by erasing the lines 
of color zero. This implies that B itself must be melonic and, due to Lemma 01 it implies that Q is unique. 
The expectation of a melonic observable B is therefore in the large N limit 



lim Ar 1 (TWT,T)) = 1 , (40) 

reproducing eq. ([TU| with 0,(13) = and R(B) = 1. Hence the melonic observables are the only observables 
of convergence order in and their expectation at leading order is 1. 

4.3.2 Arbitrary observables 

Consider now a generic observable B. The leading order contribution to eq. ([55)1 is given by the covering 
graphs Q of B having minimal degree. 

Definition 11. A covering graph of B of minimal degree Q 111111 , 

g™*\E°(g™*)=B , with uj(g min )= mm ui(G), (41) 

q ,g\e°(g)=B 

is called a minimal covering graph of B. Equivalently, the minimal covering graphs of B are the covering 
graphs having the maximal possible number of faces '^2 li F 0l {Q). 
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Thus for all B, 



lim N~ 1+Q W 

JV-s-oo 



Tr B (T,T)) = R{B) 



(42) 



with the convergence order of the observable ft(B) = ^ 2 _iy u '(G mia ) ~ W^M ^i^) an d R(J3) the number of 
minimal covering graphs of B. Take for example D = 3. Both invariants depicted in figure [7] arc of order 





Figure 7: Observables of lower order for D = 3 and minimal covering graphs. 



Q(jB) = 1 and the number of minimal covering graphs is in both cases R{B) = 3. As already mentioned, 
for matrices (D = 2) the minimal covering graphs are exactly the planar graphs with one faces of colors 12. 
Note that lemma U can be reformulated as follows: for every melonic observable B, there exists a unique 
minimal covering graph Q. 

In general determining the degree of the minimal covering graphs (hence the order Q (B) of an observable), 
and their number (hence R{B)) is a difficult problem. This is the reason for which here and below we prefer 
to treat the melonic observables and the rest of the observables separately. Indeed, in order to show that a 
tensor distributed with some /j,n converges in distribution to a Gaussian tensor wc will show that for any 
observable one can establish a large N limit 

lim N- 1+ " {B) n(Tr B (T,f)) = R(B) , (43) 

with Q,{B) and R(B) identical with those of the Gaussian distribution eq. (|42|) . However for most observable 
we will prove this indirectly, without actually computing neither Vt{B) nor R(B). It is instructive then to 
see that for melonic observables one can establish by an alternate, direct route Cl(B) = and R(B) = 1 both 
for the i.i.d. and for the properly uniformly bounded trace invariant case. 
We will need the following result. 

Lemma 6. Let Q mm a minimal covering graph of the D colored graph B with D odd (respectively even). 
Then any two lines of color ofQ nlln , l® = (v, v),^ = (w,w) <G £°(Q mnl ) share at most (respectively 
Y ) faces of colors Oi. 

Proof: Denote the number of faces of colors Oi shared by l® = (v, v) and Z° = (w, Hi) by q. We build the 
open graph Q lmn obtained from Cy mm by deleting the lines (i>, v) and (w, w) and adding four external vertices 
v, v, u> and w hooked to v, v, w and w respectively by external lines of color 0, as in figured] 

The boundary graph of Q mia , dG min is a D colored graph with four vertices. Hence it necessarily has the 
structure presented in figure [5] on the right, with q lines connecting v with w (respectively v and w ) and 
D — q lines connecting v with v (respectively w and w ). 

Consider then the graph g mnl > x obtained from Q by replacing the lines (v, v), (w, w) by two new lines of 
color zero (w, w), (w, w), like in figure [8] on the second line. It is also a covering graph of B. 

The number of faces of colors Oi of Q min and g mnl < x are respectively 

£ F 0i (g min ) = £ F°US min ) + q + 2(D - q) , £ F °i(g^x) = £ ^(^ min ) + D - q + 2q . (44) 

i i i i 

as Q min is a minimal covering graph of t/, we have 

J- p 0i (g min ) > J2 F 0l (g min - X ) => 2q < D . (45) 

i i 

□ 

Remark that this lemma also holds for D = 2. 
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w 



Figure 8: A minimal covering graph (J min ; the opened graph Q mm , the boundary graph dQ mm and Q 



Lemma 7. The convergence degree is a positive number. Moreover, for any B, there exists an infinite 
family of graphs B' such that f2(B) = Q(B'). Finally, the only normalization of the Gaussian such that both 
statements hold is the one of equation (|34[) . 

Proof: From eq. (f3U)) . Q(B) > -jyjOj(G nua ) > 0. Consider a graph B and the graph B' obtained by inserting 
a D — 1 melon (say of colors 2, 3 ... D) on one of the lines (say of color 1) of B. Call v and v the vertices of 
this melon. 

Consider a covering graph of B' , Q 1 ', such that the two vertices v and v are connected by a line of color 
in Q' . All minimal covering graphs of B' are of this kind: any covering graph of £>' such that v and v' are 
not connected by a line of color would have two lines of color sharing D — 1 faces, which is impossible 
by lemma [5] thus, 

[JJ — lj! g' ,g'\s°{g')=B' (D — 2)1 

(v,v)ee°(Q') 

By reducing the melon v and v' , B becomes B' and Q' becomes some covering graph Q of B. All covering 
graphs of B can be obtained starting from some Q' of this kind. Moreover, as B is obtained from B' by 
reducing a D — 1 dipole and Q from Q' by reducing a D dipole uj(B') = oj(B) and oj(Q') = lj{Q). Thus 

n(B) = - — ^— min u>(g)--, — 2 ——u(B) 
V ; (D - 1)! g ,0\£°(0)=8 V ; (-D -2)! v 7 

= TTT^TTj min ^0 - 777^1^0 = W • (47) 
(iJ - 1)! s' .s'\£°(s')=b' (L>-2)! 

(i>,c)e£°(e') 

By inserting D — 1 melons arbitrarily on the lines of B one then builds an infinity of graphs B' with 
Vt(B') = Cl(B). This proves the first part of the lemma. 

For the second part, suppose that one choscs a different normalization of the Gaussian measure 

( n nv dTndf " ) e ~ N " E " a T " s ^ f - ■ (48) 

n 

Then eq. (|3"5j) becomes 

7V fc(e)(D-i-^)+i- TI j4 TT T W (e) +TI j^ TT a;(e) _ ^ 
The order of convergence of an observable would then be 
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which is positive for melonic B only if v > D — 1. Moreover, if v > D — 1, then there exists only one 
observable with scaling 1/ — (D — 1), the D dipole itself. 

□ 

Different scaling of the Gaussian can make sense, but only if one decides to look at subsets of observables. 
Consider for instance a tensor with 4 indices. One can decide to only consider tensor observables in which 
the tensor effectively acts as a iV 2 x iV 2 matrix, that is the indices (1,2) and the indices (3,4) are always 
contracted between the same tensors A scaling v = 2 leads to a well defined large TV limit for these observables 
(this is just the usual large N limit of matrices). However other tensor observables do not behave well with 
this scaling: the melonic observables are arbitrarily divergent. The importance of the scaling 7V D_1 of the 
Gaussian is that it renders all the tensor observables convergent in the large N limit. 

4.4 Proof of Theorem CD 

The proof follows closely the one for matrices. Set the covariance of the atomic distribution to a 2 = 1. 
Consider the observable associated to a graph B with 2k vertices 

/i(Tr g (T,f)) = -J-^ J2 S n,nYl K -[ T ^ %-••%] ■ ( 51 ) 

n.n 7r 

Again we represent all the second order moments as dashed lines of color 0. Again we deal with the higher 
order moments in a non canonical way, by representing them as dashed lines in some pairing of T and T, 
but with further identifications one needs to track. The expectation writes as a sum over covering graphs of 
B. The trace invariant operator composes with the identifications given by the cumulants and the faces Oi 
bring each a N. One obtains 

M (Tr e (T,f)) = ]T N- k(D -^N^ F ° z ^([[d)([[K) 

g, g\£°(g)=b 

= j2 N-^-^o^^w-^o^^wqjs^qj^ (52 ) 

e, g\£ a (G)=B 

where (J\ kJ is a product over the cumulants associated to a graph. Note that if some of the lines in £°(Q) 
correspond to a higher order cumulant, all the indices of the faces Oi to which this lines belong are identified. 
We denote these extra identifications formally by the (J\ 5) . These further identifications either play no role 
(if the indices of the faces Oi on the lines are already identified in Q) , or they reduce the number of independent 
sums, hence the total scaling in N is strictly smaller than ~k(D — 1) + Y^o<i<j F^ (G) — J2o<i<j F^iQ)- 
Again V,,. , FV(G) = £ 0<i<j F^(B). 

The scaling with N of a term in this sum is therefore at most 

i - a?) + (^hr (*) < i - Jt^) , (53) 

like in eq. (|39l) . The graphs Q are covering graphs of B, and the presence of a higher order cumulant 
(potentially) brings some extra constraints on the sums. 

4.4.1 Melonic observables: direct computation 

We first consider the case when the bound ([53]) is saturated. It follows that Q is a melonic graph, and 
consequently B is melonic also. 

Two lines of color in a melonic graph Q having a unique bubble B of colors 1 . . . D cannot share all 
their D faces of colors Oi. Indeed Q is a covering graph of B, and, as it has degree zero, it is also minimal, 
thus by lemma [5] any two lines can share at most (D — l)/2 (or D/2) faces of colors Oi. It follows that all 
the lines of color of Q represent a second order cumulant (the presence of a higher order cumulant strictly 
decreases the scaling in N) 
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From lemma 01 Q is unique, thus at first order in N exactly one graph Q contributes and all lines of color 
of Q represent a second order cumulant, hence 

lim jV- 1+n(B V(Tr s (T,T) > ) = R{B) , (54) 

with n(B) = and R(B) = 1. 
4.4.2 Arbitrary observables 

We now deal with arbitrary observables. For all B only the minimal covering graphs contribute at leading 
order to cq. ([5^]) 

/i(Tr B (T,T)) = J2 N-^-^N^^^i^dl^ll + OiN- 1 )) (55) 

e min , e min \£°(e)=e 

Again, the bound in eq. (|53j) is saturated (and consequently we get a contribution from the corresponding 
gmm above) only jf the product (J\ 5) does not bring any extra identifications. Again, if Q is a minimal 
covering graph of £>, no two lines share all their faces, hence if Q possesses at least two lines coming from a 
higher order cumulant its contribution is strictly suppressed. 

Hence the graphs contributing to eq. (|55[) arc the minimal covering graphs of B such that all their lines 
of color correspond to a second order cumulant. Then 

lim N~ 1+n{B} ( Tr B (T, f)) = R(B) , (56) 

with n(B) = ■^ T y T a;(g min ) - m=2p0j(B) and R(B) the number of minimal covering graphs of B (as every 
minimal covering graph contributes exactly once), reproducing the moments of the Gaussian distribution. 

□ 



4.5 Proof of Theorem [2] 

Following the discussion of the trace invariant measures for matrices, the expectation of an observable B 
with 2k vertices for a trace invariant measure for tensors, 

Tr e (T, f )) = ]T 5 B nn £ «* [T* % k ] , (57) 



/'• 

writes as sums over doubled graphs Q D B generalizing ([25} 

' -2(D-l)fc+D|a|-ElTii p(B(a) 



t (Tr B (T,f)) 

, / p(B(a)) N 

gDB,g\£0(G)=BU^l 1 {U p ^ 1 >B p ( a )) 



N 



E(n n ^) n aiw 

ji=(«,ij)ee' [Bu^i, (uX ' e P (a)) J 

|a| 

n^(B(a),JV). (58) 

a=l 

Recall that the subgraphs with colors 1 ... D of the doubled graph Q fall in two categories. One of them, B 
(having no label a), corresponds to the initial observable, while the others B p (a) correspond to the various 
cumulants K 2 k( a ) ■ These graphs are connected by dashed lines of color and, like for random matrices, the 
Kronecker 6's compose along the faces with colors Oi. The expectation of the observable writes as a sum 
over all doubled graphs which contain B 



n 

GDB a=l 



t(Tr B (T,f)) = N- 2 ( D -V k+D W-^r( B ^)+^ FOi W fl K(B(a),N) . (59) 
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Using again the fact that the number of faces of colors Oz computes as the total number of faces minus the 
ones which don't have the color 0, the scaling with TV computes further 

H 

- 2(D - l)k + D\a\ - J2 P( B ^) + E Fij @) - E Fii{ $) > ( 6 °) 

a=\ Q<i<j 0<i<j 

taking into account that each face with colors ij, < i < j belongs cither to B or to some B p (a), 
Eo<,< 3 Fij (G) = E 0<i<j Fij @) + Elli EpS (0)) F t3 {B p (a)) the scaling computes to 

-2(1? - l)k + D\a\ - P(B(«)) + ( ° {D ~ 1] HS) + D- 

■{D-l){D-2) 1 ^~] n 1 _J2 

(I? -2; 

|a| p(b(«)) 

(D-2) 



C ( < J - '" D - 2 > *(»,(„)) + P - 1 - (Bp(tt ))) . (61) 



a— 1 p— 1 



and recalling that the doubled graph Q has 4A; vertices, = 2k while 2? has k vertices, k(E) = & and 

Elli EpS (Q)) fc(B P (a)) = fc we obtain 



M (Tr e (T,f)) = ]T n A 1^(«) 



|a| 

5DBa=l 



1 



As B and B p (a) are all the subgraphs (bubbles) of colors of the graph Q, and using lemma [1] we bound 
the scaling with N of Q by 



4.5.1 Melonic observables: direct computation 

Again we first discuss the case when the bound in eq. (|63[) is saturated. Then Q is a melonic graph such 
that every cumulant K 2 k( a ) is represented by an unique connected invariant, p(S(a)) = 1. 

As Q is melonic, B must be melonic. Furthermore Q has 4k vertices, 2k of them belonging to B and the 
other 2k to the invariants B p (a) (coming from the cumulants K2k(a))i an d ah lines of color connect some 
vertex in B with a vertex belonging to one of the B p (pi)'a. By Lcmma[5J Q is unique. Moreover its associated 
tree is the tree of B with all vertices decorated by lines of color ending in a tree vertex corresponding to 
some B p (a), hence all B p (a) = B^ 2 \ It follows that for melonic bubbles B 

lim AT-V(Tr B (T,f;f) = ( lim A"(6 (2) , N))^ , (64) 

reproducing the expectation values of melonic observables of a Gaussian distribution of covariancc K(B^) = 
limx^ K(B( 2 \N). 



4.5.2 Arbitrary observables 

Consider now an arbitrary observable B. Note that if some of the connected components B p (a) come from 
the same cumulant (p(£>(a)) > 1), the contribution of the doubled graph Q in eq. (|62[) is strictly suppressed 
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with respect to the one coming from the same doubled graph, but with all = 1. Thus at leading 

order in TV, we get 



t(Tr B (T,T) 



5DBBi(q) 



(65) 



where B\ (a) is the unique connected component of the graph representing the cumulant K2k(a) 
Among the doubled graph Q D B contributing, some represent a minimal covering graph Q 11 



decorated by a two point cumulant on all the lines of color 0. We denote such a graph Q r 



J S 0(gn 



of B 
,B( 2 ). In 



this case every cumulant is a D-dipole Si (a) = B< 2 '. We note that uj(B^) = (as the D-dipole is the first 
melonic graph with D colors). Moreover, Q mm L>£orgmin\ B^ has D 2 ~ k(Q mm ) extra faces with respect to 
gmm ^ a vj f aces f co l rs < i < j made by lines of the various B^ insertions) and 2k(Q min ) extra vertices 
(two vertices for each B^ 2 ' insertion), hence by eq. (|29p ui(Q mm Ugo/gmin) B^) = uj(Q nlln ). Separating these 
terms among the terms contributing to the expectation we get 



H Tr B (T,T) = 



E 



K{B {2 \N) 



KB) 



N 



4G n 



Rest 



(66) 



)=B 



therefore, provided that the rest of the terms are subleading in 1/N, we get 

KB) 



with 12(B) 



Um N- 1+n WjTr B (T,fj) = \ lim K(B<- 2 \N)] R(B) 

N-^-oo V / IN— too 



(67) 



lu(B) and R(B) the number of minimal covering graphs of B, reproducing 

limN^ K(B( 2 \N). 



(D-1)! W V» I {D-2)\ 

large N moments of the Gaussian distribution of covariance K(B^) 

To conclude we now prove that all the other terms contributing to the expectation are strictly suppressed 
in 1/N. Among these terms some represent non minimal covering graphs of B decorated by insertions of 
£>( 2 ' on all lines of color 0, g n - mm Ugo^n.mio) B^ 2 K For such graphs cq. holds, with w(£ mm ) replaced by 
u>(G n min ) > uj(Q nun ), hence they are suppressed. 

Consider now that Q has at least a higher order cumulant B\(a) ^ B 2 . The scaling with N of Q is, from 
cq. (IHSD (recall that p(B{a)) = 1), 



2(D-l)k+(D-l)\a\+J2F 0l (g). 



(68) 



Consider two lines of color 0, (v,a) and (a,v) touching two vertices v and v separated by a line of color j 
of Bi(a). We will compare the scaling of Q with the one of the graph Q obtained by reconnecting the two 
lines of color into a line of color 0, namely (a, a), with a B^ 2 > insertion, and reconnecting all the other lines 
touching v and v respecting the colortH (sec figure [9|). We consider the two point subgraph B^ as coming 
from a different cumulant in Q. 




Figure 9: The graphs Q and Q. 



6 If there are several lines (of colors different from 0) connecting v and v in Q, we delete them. If Bi(cc) divides in several 
connected components under this procedure, we associate a different label a to each of them (i.e. we consider each of them as 
coming from a different cumulant in Q). Both these cases give strictly subleading contributions. 
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The graph g is also a doubled graph D B, having \a\ = \a\ + 1, and F 0l (g) > F 0l {g) - (D - 1) (as 
the face of colors Oj is not affected by this change, and all the other D — 1 faces Oq touching v and v can at 
most merge two by two), thus 

- 2(D - l)k +(D- 1)H + F° l (5) < -2(D - l)k + (D- 1)\&\ + ^ F 0i (g) (69) 

i i 

and equality holds only if all the faces of colors Oq, for all q ^= j touching v and v arc merged after this 
reduction. Iterating we reduce the order of all cumulants and obtain a doubled graph representing a covering 
graph £ final of B with two point insertions B^ on all lines, £ final U £ o (gfl „ai } B^. 

At the last step we reduced a four point cumulant connected to the rest of the graph by four lines color 
namely (v,a), (a,v) and another two lines, say (b,w) and (w,b). In order for eq. (|69p to hold with an 
= sign (if not the contribution of g is strictly suppressed with respect to the one of £ final U £ o( g fi„ a i) 
it follows that the two lines of ^ final ; (a, a) and (b,b) (obtained after eliminating the insertions B^ in 
C? final U£0( g finai) B^) share all the D—l faces of colors Oq for q =^ j. Hence from lemma[H]the graph Cy final can 
not be minimal. In all cases the contribution of g is strictly suppressed with respect to the one of minimal 
covering graphs decorated by B^ insertions thus eq. (j6"T)) always holds. 

□ 
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A Perturbed Gaussian measures 

In this lengthy appendix we discuss in some detail a large class of probability measures for tensors to which 
our results apply: the perturbed Gaussian measures. This measures appear naturally in physics and describe 
random D dimensional triangulations |24[ 119) . A perturbed Gaussian measure is defined by an action 

S(T, f) = J2 T n5 Hft % + tnTrn (T, f ) 

n H 

"" = ^v)(n* D -'^)^°"' SM . PO) 

n 

with Z(tfi,N) a normalization constant. We consider only the case when all Ji are connected graphs with 
D colors, hence the most general "single trace" model. 

There are two levels of precision at which one can study the perturbed Gaussian measures. The first one 
is the perturbative level. Perturbed Gaussian measures are evaluated in terms of Feynman graphs (and the 
reader is assumed in this appendix to have some familiarity with them). The perturbative level establishes 
bounds on individual graphs. It is at this level that we work in section IA.1I . The relation between the 
Feynman graphs and the doubled graphs used so far hides a number of subtleties. The aim of the first 
part of this appendix is to first detail this relation and second to prove that the joint cumulants of such a 
distribution respect the proper uniform bounds in eq. ([7|) at the perturbative level, i.e. for each Feynman 
graph contributing to a cumulant. This result is insufficient to claim that such measures are properly 
uniformly bounded, but it constitutes a good indication that they might be. 

The second level at which one must treat the perturbed Gaussian measures is the non perturbative or 
constructive level in which one must promote the perturbative bounds to constructive bounds valid for 
the full cumulants which are sums over graphs. In order to establish proper uniform boundedness at the 
constructive level one must resum the perturbation theory. This is notoriously difficult, as the perturbation 
series is not summable (the number of graphs grows too fast) but only Borel summable (to be precise Borel - 
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Le Roy of an order fixed by the maximal degree monomial in the perturbation of the Gaussian measure). The 
resummation of the perturbative series requires a set of techniques quite different from the ones employed in 
the rest of this paper, amounting to a research field in itself: constructive field theory [55]. We will therefore 
treat in the second part of this appendix at the constructive level only a particular example of a perturbed 
Gaussian measure. The techniques we present here (generalizing the one introduced in |33| to tensors) should 
be further refined along the lines of [34) to establishing proper uniform boundedness of all the cumulants for 
any polynomially perturbed Gaussian measure. 

A.l Perturbative bounds 

We will show in the first part of this appendix that the cumulants of this measure respect in perturbation 
theory 



where the sum runs over all -D-colored graphs B with 2k vertices and p{B) connected components B Pl and 
lim.Ar-i.oo K(B 1 N) exists and is finite for all B. By theorem [2] all the perturbed Gaussian measures become 
Gaussian in the large N limit in the perturbative sense. It is however naive to conclude that the large N 
limit of such models is trivial. The covariance of the large N Gaussian (equaling the large N expectation of 
the D-dipole observable B^) is the full resumed two point function of the model, and has a very non trivial 
dependence on the parameters tu |20j . 

We will evaluate the moments and cumulants of the measure (|70p by expanding in Taylor series with 
respect to t-u (hence the name "perturbative results"). As we will see in the second part of this appendix 
the ensuing series are not summable, and in order to promote these results to the non perturbative level one 
needs to work quite a lot. The joint moments of the probability distribution of tensor entries 



are expressed as sums over Feynman graphs. They are obtained as follows: upon expanding with respect 
to tn, each invariant Tr%(T, T) (represented by a graph H with D colors 1,2... D) becomes an insertion in 
a Gaussian integral. The Gaussian integral is then evaluated in terms of contractions (pairings) of tensor 
entries. For each such contraction scheme one draws a Feynman graph. The invariants act as effective 
vertices (interactions) of the Feynman graphs (not to be confused with the black and white vertices of H 
itself which represent the tensor entries T and T). The effective interactions are supplemented by effective 
lines, (propagators, Wick contractions) representing the pairing of two tensors and with the Gaussian 
measure. We represent the contraction of two tensor as a dashed lines of color connecting the corresponding 
black and white vertices. Thus a Feynman graph Q has D + 1 colors, for the dashed lines and 1 . . . D for 
the effective interactions. 

The insertions , T= i . . . T=* in the joint moment are represented as external black or white vertices 
of valence 1. The external vertices are joined by lines of color to the rest of the Feynman graph. Thus 
the dashed lines of color fall into two categories: internal joining two tensors (that is black and white 
vertices) on two effective interactions % and %' and external joining an external vertex with an internal 
vertex on some H.. The Feynman graphs are then are nothing but open D + 1 colored graphs of definition 
IH1 Two examples of Feynman graphs are presented on the left in figure The effective interactions H are 
represented with solid lines of colors 1,2 and 3. Both graphs have four external lines of color 0. 

The cumulants «2fc (T^ , . . . T» J are sums over connected Feynman graphs Q with 2k external (uni- 
valent) vertices. We stress that the G's contributing to a cumulant are connected as a graph with D + 1 
colors. 

Each D + 1 colored graph represents an abstract D dimensional simplicial pseudo manifold [17] . This 
pseudo manifold is obtained by associating a D-simplcx to each (black and white) vertex in the graph (hence 
to each tensor entry T and T). The D — 1 simplices bounding the D simplex are colored 0, 1 up to D. 
This induces colorings on all lower dimensional simplices. The D simplices are then glued respecting all 




fc(B)=, 



(71) 




(72) 
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the colorings: a line in the graph represents the unique gluing of two D simplices along boundary D — 1 
simplices which respects the colorings of the D — 1, D — 2 etc. simplices. An effective operator Tr^(T, T) 
with 2k tensor represents the gluing of 2k simplices around a vertex forming a "chunk" . For example in 
three dimensions an operator represents a gluing of tetrahedra around a vertex. The boundary of such a 
chunk is paved by triangles (represented by the half lines of color 0). Topologically the chunks are cones 
over their boundary, hence they can have non trivial topology. A Fcynman graph represents the gluing of 
such chunks into a pseudo manifold. As the combinatorial weights and amplitudes of the graphs are fixed by 
the Feynman rules, the measures (|70p encode a canonical theory of random pseudo manifolds in arbitrary 
dimensions. The leading order melonic graphs represent spheres. 

Each contraction in the Gaussian integral (hence dashed line of color 0) replaces the two tensors and 

by a covariance N h-i &np- The amplitude of a Feynman graph is then 

ag = E(n^)^- x e (e) )( n ^rk, 

n,n H(G) l°=(v,v)eS°(g) i=l 

u{g) 

E(n n *«*«0( n nM< <™ 

where T-L{G) runs over all the subgraphs with colors 1 . . . D of Q, H(Q) = \H(G)\ denotes the number of 
such subgraphs and |£°(£?)| the number of lines of color of Q. The Kronecker S's compose along the faces 
with colors Oi. The faces with colors Oi of Q are either internal or external (see section B~2"]) . The internal 
faces Jfat(&) ( w ith ^iSt(^) = l^hit^(^)l) yield a f ree sum hence bring a large factor N. The external faces 
/ G -^ext necessarily start and end on two external vertices u and u corresponding to two external insertions 
T and T in the joint moment, / = [u, u). The amplitude of a graph becomes 

ag = (n^( e ))^ iweK( ^ i)r(e)i+E *^ (e) n *nt«i, (74) 

•HG?) /=(«,s)6u,^ t (e) 

and the operator 0/=(u crteUj.?" * rg) S n i a » reproduces the trance invariant operator associated to the bound- 
ary graph dQ (see again section H~2")) . We finally conclude that the amplitude of a graph is 

As already mentioned (the second example in figure [5]) in spite of the fact that Q itself is connected, the 
boundary graph dQ can be disconnected. It follows that a cumulant, which is a sum over connected graphs 
Q expands as a sum over all possible D colored graphs (connected or not) corresponding to the possible 
boundary graphs B = dQ 

aCT^f* ...t^) = £ [ £ JL ([] t mg) ) ^W-M^I+a^W)] , ( 76 ) 

where B runs over all closed D colored graphs with 2k vertices (connected or not), and Q runs over all D + 1 
colored connected graphs with 2k external legs of color whose boundary is B. For every Q, T-L(G) runs over 
its internal subgraphs with colors 1, 2 up to D, \£°(G)\, Ff£ t (G) and H(Q) denote the total number of lines 
of color (including the 2k external lines), internal faces of colors Oi of G, and respectively subgraphs T~L. 
Finally, S(G) is some symmetry factor. 

We can now describe the precise relationship between the Fcynman graphs and the doubled graphs used 
to establish theorem [2] The doubled graphs for a perturbed Gaussian measure consist in the observable B 
and the boundary graphs B p (a), p — 1, . . . p(B(a)) of the various Feynman graphs G{o) contributing to each 
of the cumulants K 2 k( a ) arising in an expansion in cumulants of the moment ^(Trg(T, T)). 
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In order to conclude at the perturbative level that all perturbed Gaussian measures become Gaussian in 
the large N limit, we must prove that for all Q with dQ = B = U^ff^Sp 

i y(r>-l)ff(S)-(D-l)|£ (S)|+Ei^ »i(e) < N -2(D-l)k{dG)+D~p(dG) ^ fnjx 

that is the cumulants respect eq. (|7Tj) . 

Theorem 4. For every connected D + 1 colored graph Q with 2k external vertices, \£°(Q)\ lines of color 0, 
-Fint(^) internal faces of colors Oi, H(Q) subgraphs with colors 1,...D and p{dQ) connected components of 
the boundary graph dQ 

(D - l)H{Q) -(D- l)\e°{9)\ +J2 F ^(G) < -2(X> - l)k(dG) +D- p{dQ) . (78) 

Proof. The proof of this theorem is divided into two parts. We first present an iterative algorithm which 
reduces the graph Q to the D + 1 colored graph dQ U £® xt (G) consisting in the D colored graph dQ decorated 
by an external line of color for each of its 2k vertices. At each step of this algorithm we will obtain a graph 
Q("> interpolating between Q^> = Q and Q( s ™^ = dQ U £® xt (Q). Second we will prove that at each step of 
this algorithm the quantity 

Q(s) =D- C(0«) + (D- l)[H(Q^) - C(gW)] - (D — + J2 F "nt(G {s) ) , (79) 

i 

with C(Q^ S ') the numbers of connected components of Q^ s \ H(Q^) the number of bubbles (subgraphs) with 
colors 1, 2 up to D of Q^, \£ °(£ (s ))| the number of lines of color of Q^ and i^ t (S (s) ) the number of 
internal faces of colors Oi of^( s ) is strictly increasing. As Q(0) = (D - l)H{Q) -(D- 1)|£°(<7)| + X)» F int(G) 
and Q(s max ) = D — p(dQ) — 2{D — 1) k(dQ), we conclude. 

Obtaining dQ U £e Xt (£). The algorithm we present here has been introduced in [32] . 

Consider a connected D + l colored graph Q^ with 2k external vertices. We first define an internal q+1 
dipole with colors 0, ii, . . . i q as two internal vertices v and v of Q^ connected by an internal line of color 
and exactly q lines of colors ii, . . . i q . An example of an internal q + 1 dipole with colors 0, 1, ... q is given 
on the left in figure 1101 An internal q + 1 dipole can be contracted. The contraction consist in deleting the 
two vertices v and v and the q+1 lines connecting them, and reconnecting the remaining lines respecting 
the colors. 




Figure 10: A q+1 dipole with colors 0, 1, ... q. 



Under a contraction we obtain a new graph C/( s+1 ) having two less vertices, one less internal line of color 
0, q less internal faces of colors Oi and the same number of external vertices, 2k. Note that the new graph, 
g(s+i)^ can potentially be disconnected. Note also that neither the external vertices of Q^ s \ nor its internal 
vertices hooked by a line of color to external vertices can be deleted. 

Consider the graph obtained starting from Q and contracting iteratively, in an arbitrary order, the 
maximal number of internal q + 1 dipoles with colors 0, i\, . . . i q . The number of internal dipoles contracted 
equals the number of internal lines of color of Q, \£i nt {Q)\- We show below that the final graph 5' Sm "' 
is dQ U £%xt(G), the boundary graph of Q decorated by an external line of color on each of its vertices. 
Examples of this full reduction are given in figure 1111 

The final graph Q( s ™^) has 4k vertices, 2k coinciding with the external vertices of Q and 2k with the 
internal vertices of Q hooked to external vertices by external lines of color 0. It has no more internal lines 
of color but still has 2k external lines of color 0. As the internal vertices are each touched by exactly one 
line for every color 1, 2 up to D, Q( s ™»^ has exactly k lines of every color 1, 2 up to D. Furthermore Q( s ^^> 
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Figure 11: The reduction of all the internal dipoles in a graph. 



has no internal faces of colors Oi. However the external faces with colors Oi can never be deleted by this 
procedure hence all the faces of colors Oi of C?( Smax ) are external and they are one to one to the Dk external 
faces of colors Oi of Q . It follows that all (external) faces Oi of (y( s >»>»<) contain exactly one line of color i, 
connecting the two internal vertices hooked to the external vertices which share the face Oi. The lines of 
color of g( s »-) are all external lines and one to one to the external lines of G, S°{Q l - s ' a ^') = £® xt (G). By 
deleting the lines of color (and flipping the black and white vertices), the final graph Q( s ™-^) \ £°(^G^ Sm ^) 
will have a vertex for every external point of Q, and a line of color i connecting two vertices u and u for 
every external face / = (u,u) of colors Oi of Q. Hence £/( s ™"<) \ £°(C/( s ™ e "')) = dQ. 

Bounds. Suppose we reduce a dipole of colors 0, 1 . . . , g to pass from (/ W to C/( s+1 ) . We have two cases. 
Either the two vertices v and v belong to two different bubbles (connected components) with colors 1, 2 up 
to D and the dipole is necessarily a 1 dipole made exclusively by a line of color 0, or the two vertices belong 
to the same bubble with colors 1. 2 up to D. 

First case. We have v G Hi and v G %2, and both Hi and H2 belong to the same connected component of 
G^ . The number of connected components does not change by contracting the dipole, C(Q( S+1 >) = C(Q^ S '), 
To see this, consider the bubble Hi- As it is a graph with D colors it can not become disconnected by deleting 
v. Chose a spanning tree T-y in Hi \ v (the bubble with v omitted), and a tree T2 in H2 \ v. Complete it 
by adding the lines of color 1 touching v G l\ and v G l\ and the line of color = (v, v), and finally to a 
tree in the entire connected component of G^ by adding lines T re8t . The tree T\ U l\ U U^UT 2 U r rcst 
becomes after reduction the tree T\ U l l UT2 UT rest (with I 1 the new line of color 1), spanning one connected 
component in Q( s+1 \ 

The two bubbles Hi, Hi C Q (s) are collapsed into a unique bubble of£< s+1 ) thus H(g ( - S + 1 ^) = H(G^)-1. 
The number of lines of color decreases by 1, \£ °(Q( S+1 >)\ = \£°(G^)\ — 1, and the number of internal faces 
of color Oi docs not change ^°U5 (s+1) ) = ^°U^ (s+1) ) hence 

D - C(£ (s+1) ) + (D - l)[H{g( s+1 '>) - CiG^)] - (D - l)|£°(g (s+1) )| + F?nt(G {s+1) ) 

i 

= D- C(g&) + (D- l)[H(G^) - C(G^)] -{D- 1)\£\G {S) )\ +^if n i@ (s) ) ■ (80) 

i 

Second case. Both v and v belong to the same bubble v,v G H. In this case the number of connected 
components of Q( s > can increase when contracting the q + 1-dipole (note that, like in the previous case q 
can be zero, but it can also be larger than in this case). As each of the new lines D — q lines (one for 
each color not belonging to the q + 1 dipole) must belong to some connected component of G^ s+1 \ we have 
C(g( s+1 > ) — C(G^ ) < D — q—1. Moreover, if one of these lines belongs to a connected component created by 
the contraction, then it certainly belongs to a new bubble of colors 1, 2 up to D created by this contraction. 
Hence C(G {s+1) ) ~ C(G {s) ) < H(G {s+1) ) - H(G {s) ). As before, \£°(G {s+1) )\ = \£°{G {s) )\ - 1, but q internal 
faces of colors Oi are deleted, ^ F^ t (G {s+1) ) = F ?nt(G {s) ) ~ 9; hence 

D - C(G [S+1) ) + (D- 1)[H{G {S+1) ) - C{G {S+1) )] -(D- l)\£°(G {s+1) )\ +Y, F ?nt(G {s+1) ) 

> D-C(G (s) )-(D-q-l) 
+(D-l)[H(G {s) )-C(G {s) )] 
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-{D l)\E°{g^)\ +D-1+Y, & {S) ) <? > (81) 

i 

thus in both cases Q(s + 1) > Q(s). □ 
A. 2 Prom perturbative to constructive bounds 

In order to prove that the full resumed cumulant is properly uniformly bounded one must perform a good 
deal of extra work. Before proceeding to the non perturbative treatment, we first establish a technical lemma. 

Lemma 8. Let a and r and £ be three permutations of k elements. We denote c(r) the number of cycles of 
the permutation r. Then 

c(0 + c(a£) + c(r) + c(gt- x ) < 2 + 2k . (82) 

Proof: To triple of permutations £, a and r we associate a ribbon graph constructed as follows. 

Consider the decomposition of £ in cycles, £ = C\ . . . C c m, each of length \C r \. We draw a fat vertex for 
every cycle of £ having A\C r \ hairlines. We assign a label l q a q f3 q j q , ^(g)a{(g)/%ij).7£(ij)> ■ • ■ turning clockwise, 
to each half line, see figure[12] Two consecutive halflines share a strand. The strands are of three kinds: solid 
(connecting j q to l^t q )), dashed (connecting a q with (3 q ) and wiggly (connecting l q to a q or j3 q to j q ). Thus the 
halflines representing Vs and j's are solid -wiggly, and the ones representing a's and /3's are wiggly-dashed. 
We represent the permutations a and r by ribbon lines connecting the halflines l q to j a ( q ) and ct q with /3 T ( g ). 
The ribbon lines l q — > j a ( q ) arc then solid-wiggly and the ribbon lines a q —> P T ( q ) are wiggly-dashed. 




Figure 12: Ribbon graph associated to three permutations £, a and r. 



The faces (closed strands) of the ribbon graph thus obtained are of three types: 

• dashed faces (5 q — > a q — >• P T (q) ■ ■ ■ ■ We call them "r faces" as they are indexed by the permutation r. 
Their number is the number of cycles c(r). 

• solid faces j q — > l^ — > j^( q ) We call them "tr^" as they are indexed by the permutation er£. 

Their number is the number of cycles c(cr£). 

• wiggly faces f3 q — > a T -i( g ) — >■ l T - 1 (q) ~^ iaT- 1 {q) ■ ■ ■ ■ We call them "<tt _1 " faces as they are indexed by 
the permutation ot^ 1 . Their number is the number of cycles c(ar^ 1 ). 

The ribbon graph has c(£) vertices and 2k lines, hence 

c(£) - 2k + c«) + c(r) + c{ot- x ) < 2 . (83) 

□ 

We are now going to the core of the non perturbative definition of the simplest quartically perturbed 
Gaussian measure. We will denote the indices either as n or as na with n = n 1 and a = n 2 , . . . n D . We will 
deal with the measure 

S(T, f) = 2 ^SnnS^nS + A E T - f f-'V fil 5 nn<5 a a<5pp^| 
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^=^(n» D -«)-™°-' s< ">. <w 

flQ 

The quartic perturbation corresponds to a melonic invariant. From now on all the repeated indices are 
summed. We consider the generating function of the moments of our distribution 

Z(J,J-,X,N) = J(l[N^ d ^) e-^^^-W^) . (85) 

na 

For every finite N, Z(J, J; A, N) is denned by an integral over N D complex variables which is absolutely 

convergent (and bounded by e ^ n & for 3? A > and divergent for 5f A < 0. Thus Z(J, J; A, TV) is 

an analytic function in the right half complex plane. The generating function of the cumulants (connected 
moments) is W(J, J; A, N) = In Z(J, J; A, N), 

^n A ,% lSl ,...T nkSk ,f nkiSk ) =N -mn-i) & 2k) W(J,J;X,N) (86) 

UJ n 1 3 1 UJ fil g 1 . . . UJ nkak U •Jfi k a k j-.j—v 

The perturbative treatment described in the previous section consists in two steps. To compute the 
cumulants (connected moments) one first takes the logarithm of Z( J, J; A, N) and second one expands this 
logarithm in Taylor series around the point A = 0. Both steps are in fact problematic. While one has a 
well controlled expression for Z{ J, J; A, N) as an absolutely convergent integral, the same does not hold for 
W(J, J; X,N). Furthermore, A = belongs to the boundary of the analyticity domain of Z{J, J; X,N). So 
far we don't know anything about the analyticity domain of W(J, J; A, N), but A = can certainly not be 
in its interior: it can at best belong to its boundary. A Taylor expansion around a point belonging to the 
boundary of the analyticity domain of some function typically leads to series which are not summable, but 
only Borel summable. 

One needs to find an intermediate expansion of W( J, J; A, N) (hence a proper way to take the logarithm) 
which yields a definition of W(J, J; A, N) in some analyticity domain, such that A = belongs to its boundary 
and furthermore W(J, J; A, N) is the Borel sum of its Taylor expansion around A = 0. Moreover, as we are 
interested in the N — > oo limit, both the convergence of the constructive expansion in its analyticity domain 
and the Borel summability around A = must be uniform in N. 

This intermediate expansion is called the constructive expansion. The constructive expansion captures 
some features of the perturbative expansion (for instance the perturbative bounds on Feynman graphs can 
be promoted to bounds on the terms in the constructive expansion) but unlike the former it is absolutely 
convergent and provides the good non perturbative definition of the connected moments of perturbed 
Gaussian measures. 

The draw back of the constructive expansion is that it is rather involved and its study requires a number 
of new techniques. Note that a priori one can give several constructive expansions of the same measure. 
In particular the constructive Loop Vertex Expansion (LVE) introduced in [33] has proved (the only one) 
adapted to matrix models. We will extend in the sequel the LVE to the random tensor model with measure 
given by eq. 



A. 2.1 The Loop Vertex Expansion 

The LVE is a combination of three ingredients: the Hubbard Stratonovich intermediate field representa- 
tion, the universal Brydgcs-Kcnncdy-Abdessclam-Rivasscau forest formula and a replica trick. The last two 
ingredients are a recurrent feature of any constructive expansion while the first one is specific to the LVE. 

Step 1: We use the Hubbard Stratonovich intermediate field representation of the interaction term 

ab 
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where a a b is an auxiliary TV x N matrix field. Denoting I the identity matrix of size N D 1 x N D 1 we write 
more compactly 

ab 

thus Z( J, J; A, N) becomes, denoting 1 the identity matrix of size N x N. 

z «j;A, W )=/(n^-«)(n^) 

na ab 

-<?ab<fab-N D - 1 f.(im + J- 1 ^ T (<T-<T'<)<®l) T p3 +JV D_1 f |J J+JV"" 1 ! jj j 

e nP \ v N ' np iP a np nl3 . (89) 



As o~ — <7^ is antihcrmitian, 1 + ■>/ N u-i (& — o^) is invcrtiblc. We call its inverse the resolvent R(o~) 



-l 



l + i/ n d-i (c — ^ ) and note the following properties 



-trln(i?(a)) = -J j^R pn , — ptrln^)) = J j^Rpn , 



\\R(a)\\<l VA>0, (90) 

where ||i?(cr)|| denotes the operator norm. The introduction of the intermediate field renders the integration 
over TT Gaussian, thus 

Z(JJ-XN) = [ (TJ d(Tfl P ddfl P j c ' trCTCTt + Trln ( fl ( g )^ I )+^ P ' 1J n fl( fl ( CT )'» 1I )^«. fJ e J pg 

J np ^ 7Fi 

_ f (Y[ ^ g "P^"P ^ c -trg-g- t +jV c ~ 1 trln(fl(g))+jV p ~ 1 tr(fl(a-)J,/ t ) ^g-^ 
J np ^ m 

where tr denotes the trace over an index of size N, Tr denotes a trace over an index of size N D , and 
(J J')p n = JpBib~g a jna is a TV x N hermitian matrix of external sources. Note that J and ,P are independent. 

Step 2: The second ingredient consists in evaluating the integral over a by a replica trick. Let X be a complex 
vector of components X\ , . . . ■ We want to compute an integral with normalized Gaussian measure of 
covariance C (denoted dfj,c(X)) of some perturbation V(X, X). We expand in V (of course all this is justified 
only provided that the final expression is absolutely convergent) to get 

/ = f dfi C (X) e v ^ = f df,c(X) £ 1 V(X,X) n ■ (92) 

J J ri>0 n ' 

The term of degree n can be rewritten as a Gaussian integral over n replicas X^\ . . . X^ with degenerate 
covariance between replicas C^V^ = C a i, hence 

I = Y J -,jd^{X^ l \...X^) Y[V(X^,X^). (93) 



n>0 



We will regard each term in this expansion as a function of parameters w tJ , evaluated for all w lJ = 1, 
corresponding to a Gaussian measure with covariance C^g = C a i , C^j = w l ^C a i j. 
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Step 3: The third ingredient is the universal Brydges-Kennedy-Abdesselam-Rivasseau forest formula [3"5] . 
Consider n sites labeled 1, 2 ... n and a function / depending on link variables i y with i ^ j . Then 



,a,..i)-£jf( n *«)(n^ 



W kl (F n ,u) = 



inf 



(94) 



where J-„ runs over all the forests (i.e. graphs with no loop lines) with vertices labeled 1, 2, ... n built over 
the n sites, |£(.F„)| denotes the number of lines in the forest T n , Vk^iiFn) is the unique path in the forest 
T n joining the vertices k and I, and the infimum is set to zero if there is no such path (i.e. k and I belong 
to different trees in the forest). 
We compute 



n 



£>ab;i,j x ' 



dx l i 



d^ JC jx^\...x^) Y[v(x^,x^) 

1=1 



-c 



n (e 



6 c- 6 



?8Xi^ ab 5X^ 



n 



1=1 



[d^ c jxu,...xM)( n (T,^m c ^^w))f[nx^,x^ 



(95) 



(t,3')e.F» ah """" """-6 

hence the forest formula applied to the replicated integral yields (the repeated indices a and b are summed) 



'=e^e A n duij ) fd^ iM cjx^\...x^) 

n 

Y[v{x^,x^) , 



6 c- 6 



(96) 



with w u (J r n , u) = 1 and (J-'n, u) = m£/).,l)EVi-> (:f„) ufe/ - Thus in our case we get (taking into account that 
the measure over er is 1(8)1) 

Z(J,J;X,N) = EAEA II [d^HFn,u)im(<r {1) ,---° (n) ) 

n 

Y[{N D -hrln[R(a^)] + N D -Hr[R(a^) JJ+] } . (97) 



n 



Note that the Gaussian integral factors over the connected components of the forests (i.e. trees). The 
main advantage of eq. (|97| is that it allows to compute W( J, J; A, N) very easily: whenever a function is a 
sum over forests of contributions which factor over the trees, its logarithm is the sum over trees of the tree 
contribution, hence 

W(J,J;X,N)=J2^] J2 [ ( II dulJ ) f d ^(T n ,u)im(<r (1 \---<r (n) ) 



n>l T, 
S 



(ij)er„ ° a ab ba - 



(98) 
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where T n runs over all trees with vertices labeled 1, 2, ... n and w" (T ni u) = 1 and w l ° ' (T n , tt) = mi^^^'p (j- n 
with Vi^,j(Tn) the path in the tree T n connecting i and j. Expanding the product, we get 

W(J, J; A, N) = 

5 S \ 



,kl 



(99) 



(1) 



) n 



EII[ JJt W d E ^ fc) u^ 



2h 



fc=0 d=l 



«l<«2-"<ifc 



[] trln^^W)] 
5 5 



i = l 



=E^ (D ~ i)n E/ 1 ( n **) l*wK,«)v*x{o*\...oM)( n 

«>i ' r„ J ° (i,j)er n J V(i,j)e- 

ti 1 ft; n n 

EryIl[ Jjt ]^n d E % (il) k fll % (i2) U-^ (ii) )^ II trln[ii(a«)] . 



•l,<2---.«fB = 1 



i=l 
i^i!,...i k 



We represent every vertex of 7^ corresponding to a trln[i?(crW)l as a fat vertex, and every vertex corre- 
sponding to a term R(cr( ir ') 

n r m r 8j flattened vertex (whose end points we label n r and m r ). We designate 
the flattened vertices as external. The vertices are labeled by the index i of the corresponding replicated 
field o-M. 

Each derivative with a and a' brings a resolvent (of the appropriate replica). The resolvents are contracted 
along the lines of the tree which (as the field a has two indices) are double (ribbon) lines. The two indices 
of a provide a well defined ordering of the tree lines touching each vertex generating, for each 7^, a sum 
over the various unrooted labeled plane trees T^ k with n vertices, out of which k (having labels ii,...ik) 
are external, compatible with T n . In figure [13] we presented such a plane tree having two external vertices, 
1 and 4 (we added a dashed strand for the flattened vertices). 




Figure 13: A labeled plane tree in the LVE. 



Every tree T n with assigned degrees of the vertices di,...d n , has ! . . . di k ! Ilt^* (^« ~ 1)' associated 
plane trees 1~Pu , corresponding to the permutations of all but one of the halflines touching each vertex (and a 
choice di r of where to place the external insertion on the external vertices) . As the number of combinatorial 
trees with assigned degrees d\ , . . . d„ is ^ dl _iy~ 2 ^ Ziyi we get 



E*= E 



(n-2)! 



d 1 .-.-d n = l 

J2 di=1n-2 



(d x -l)\...(d n -l) 



——>d,.\...d,\\\ul, 1: 



= (n-2)! J2 



d ik = (n-2)\ 



d 1> ...d rl = l 
Y,di=2n-2 



2n + k-3 
n-2 



as the sums over di yield the coefficient of the term of degree x 2n 2 in the expansion of 



1 



1-x' 



l — k 



(1-x) 



„\n—k 



(1-x) 



Tl+fe 



= * n E 



n + k + p — 1 



(100) 



(101) 
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The sides of the ribbon lines and of the fat vertices form open strands. The tree in figure [T3] has two 
strands. To every plane tree with external vertices one canonically associates the ordered lists of vertices 
encountered when following the strands starting from their external index n q of the vertex i q (including the 
two end vertices of the strand). For the example of figure [T31 this lists are 1, 2, 3, 2, 4 and 4, 2, 1,5, 6, 5, 1. A 
moment's reflection reveals that the contribution of each plane tree is proportional to the ordered product 
of resolvents associated to the vertices along the strands: for the example in figure [T3] it reads 



R(aW)R(aW)R(pW)R(o-W) R (o-W) 



r(2 



(3> 



r (2)> 



r (4)> 



R(a^)R(a^)R(aW)R(a^)R(a^)R(a^)R(a^) 



r(2h 



r(D 



(5)> 



r (6)> 



r (5)^ 



We index the strands of T^ k by their start and end vertices i& and i(_ 



n 2 ni 

We denote Eld 
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4 £(d)- vve uexiuue ^[u — t £,(d)] the 
ordered list of of the vertices encountered following the strand (including the external vertices) starting from 
the index m d of id and ending on fh^u) on As the strands are the boundary of the plane tree, £ is a 

cyclic permutation (each cycle of £ corresponds to a connected component of 7^.). Derivatives w.r.t a and 

<jt bring factors — \ I n d-i an d \ n d-i respectively, and we finally obtain 



W(J, J; A, N) 



N 



D— 



n n .1 

n>i k=o n,i 2 ...,i fc =i r o JU (ij)er„ 



/ 



jes[d->£(d)] 



d=l 



A r 



D— 



'Ea^Es E E/( n 

«>i fc=o i 1 ,i 2 ...,i fc =i r o (jj) € r„ 



du 13 



d=i 



n ^ 

jGH[5<i-i(l)^ d (l)] 



(103) 



where T n is the unique combinatorial tree to which the plane tree T^ k reduces. The product over d is the 
ordered product of resolvents and external sources encountered when going around the tree. We chose as 
start point of this product the vertex i\ but, as the trace is cyclic, one can chose any other vertex 12, 13 and 
so on as the start vertex. We can now prove our first result concerning the LVE expansion. 

Theorem 5. The series in eq. (|103[) is absolutely convergent for < A < 3~ 2 and \\JJ'\\ < 3 _1 (hence 
uniformly in N). 

Proof: We bound tr(rj^i) < ivrjll^ill, and take into account ||iZ(cr)|| < 1 for A > 0. The Gaussian 
integrals arc normalized, and the integrals over the parameters u arc bounded by 1, thus 



n 

\W(J,J;\,N)\ < ^E^ A "^EM l|Jjt|1 



n>l 



n 

E 

> lt i 2 ...,i fc =i 7-0 



Ei 



= Y- 



k=0 



\\JJ f \\' 



2>1 



k=0 



n\ , (2n + k- 

Jn—k)\ {n - 2) \ n-2 

(2n + fc-3)! 
k\(n - k)\(n + k - 1)\ ' 
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and using -r 



(2n+fc-3)! o2n+fc-l t 

!(n-fe)!(n+fc-l)! < d W(3 § et 



1 ' . 

\w(j,j;\,n)\ < ^ D E( 32A rE( 3 ii JJt n) . 



(105) 
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which converges for < A < 3 2 and 1 1 J,P 1 1 < 3 1 . 

□ 

Equation (|103[) yields the non perturbative definition of the generating function of the cumulants of our 
probability measure (in the interval A € [0, 3 -2 ) of the real axis). In order to show that W( J, J; A, N) thus 
defined is the Borel sum of connected Feynman graphs (i.e. it is the Borel sum of its Taylor expansion around 
A = 0) we must show that (|103|) is a Borel summable function in A uniformly in N. 

Theorem 6 (Nevanlinna-Sokal). A function f(X,N) is said to be Borel summable in A uniformly in N if 



f(X,N) is analytic in a disk !JJA > R with R independent of N 



• /(A, N) admits a Taylor expansion at the origin 

r-l 

f(X,N)=Y / fN^ k + R N ,rW , I^W(A)| < Ko- r r\\X\ r (106) 

k=0 

for some constants K and a independent of N . 

If /(A, N) is Borel summable in X uniformly in N then B(t, N) = J2T=o h.f N M * s an analytic function 
for \t\ < a" 1 which admits an analytic continuation in the strip {z\ \§fz\ < cr -1 } such that \B(t,N)\ < Be 1 !** 1 
for some constant B independent of N and f(X,N) is represented by the absolutely convergent integral 

1 f°° 

/(A, N) = - jf dtB(t,N)e-* (107) 

That is the Taylor expansion of /(A, N) at the origin is Borel summable, and /(A, N) is its Borel sum. 

Theorem 7. The function N~ D W{J, J; A, N) with W(J, J; A, N) defined in eq. (|103l) is Borel summable in 
X uniformly in N for 1 1 J J' 1 1 small enough. 

Proof: Consider a complex A = |A|e lv . Following step by step the proof of theorem[S]and using ||l + /oe M || > 
|sina| =*> \\R(o-)\\ < 7— r it follows that (fTU5)| is convergent for |A| < 3 _2 | cos f | 2 = costp+1 , hence it 

|cos(f J| 

certainly converges in the Borel disk 18 < 5RA -1 . 

To compute the remainder i?Ar r (A) we separate N~ D W( J, J; A, N) into two terms: the terms with 
n < r + 1 and the ones with n > r + 1. The terms with n > r + 1 arc all in the reminder and admit the 
bound 



n=r+l k=0 i 1 ,i 2 ...,i fc =l T JU (i,j)eT„ 

k 

/^(r,„u)i^i^ (1) ,-.^ (n) )tr[n(jjt 11 fl^Ci))) 



du l ° (108) 



jesK"-i(i)->e <i (i)] 



< E iai^Eii^ii* kHn ml+k-iY - S lAr'Eii^ii^+^iAr^. 

n=r+l k=0 ' V n=r+l k=0 

with K some constant and both ||JJ^|| and A small enough, which is certainly bounded by _ftT r r!|A| r . It 
remains to find a good bound for the contribution to the reminder of the terms with n < r + 1 , 



T' 71 71 2. 

^e^-^Ef e e/( n 

n=l ' fc=0 ' ii,<2-,» fc =lr° t (i,i)£T„ 

y^^,^!^ 1 ),...^))^^^ n • (no) 



d=l ieHK d-l (1) ^d (1)] 
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For each plane tree, we use a Taylor expansion in with integral reminder of the product of resolvents up to 
some order to be chosen later 



(in) 



The first terms yield some series in A, as the Gaussian integral is non zero only for an even number of 
insertions. For every resolvent appearing in the product we have, taking into account that a — commutes 
with i?(cr), 



1 



■^^^[R(ta)(a-^)R(ta)} ab =(a np ^ + at p -^)R(ta) ab , (112) 



where the (repeated) indices n and p are summed. Taking into account the copies we get, 



o t tr[n(y n R ^ uy 



r Wt. 



up (i) »np (i)t 



tr 



[n 



jesK"-i(i)^"(i)] 
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Integrating by parts the Gaussian integral we get 



( n 

(i,j')er„ 



(T„,u)l®l 



(aW,...a(»))tr[ f[ W f ( II 



f ( II duy) [d» wij(TntU)im (aW,...aW) (2V 



dec(5) jes[d->-e(<Q] 
S S 



n % 3. 



W)er„ 



<5ct W <5<7 0)t 



tr [n( jjt n 

d=i 3 - eH [^-i(i)^-e d (i)] 



(114) 



It follows that the derivative of order s is 

r-l 



QS 

dt s 







n ^ 

(«,i)er„ 



/^« ( r B ,„)i®i( CT (1) J ..^ (n) )tr[ll(jJt n i?(t CT «; 

d=i jesK"-! (1)^(1)] 



5 \« 



('j)er, 



tr [n( jjt n 



(115) 



When computing explicitly the derivative operators acting on the trace one generates ribbon loop lines 
decorating the plane tree T n .k- The traces recomposc to reconstitute the product of R(cr) and J J' on each 
face of this graph. An example is presented in figure [14] consisting in the tree of figure [13] decorated by two 
loop lines. Its contribution is 



R(taW)R(taW) tr R(ta^)R{ta^ 
JJ^R(ta^)R(ta^)R(ta {3) )R(ta^)R(ta {4) )JJ^ 

i?(ta( 4) )i?(ta (2) )i?(ta( 1) )i?(t < T( 5 ))i?(i ( 7( 6 ))i?(t ( 7( 5 ))i?(i ( T( 1 )) 



(116) 



36 



Figure 14: A plane tree decorated by loop lines. 



Bounding again the resolvents by 1, and the traces by N times the norm, each such term is bounded by 

^) S ||Jjt|| fc A^, (117) 

as the number of faces of the ribbon graph obtained from the plane tree T^k by adding the s loop lines is at 
most 1 + s. 

The initial tree has 2(n — 1) + k resolvents. Every derivative brings a new resolvent, hence the number 
of contractions (the number of ways one can connect the loop lines on the tree) is 

[2(n-l) + fe][2(n-l) + fc + l]...[2(n-l) + fc + 2s-l]= 1 p(n - 1) + fc - 1]! ' (118) 

Choosing s = r — (n — 1), and taking into account that w lJ < 1, the Gaussian integrals are normalized and 
the integral over dt is bounded by 1, the remainder term is bounded by 



"-'£>r'X;^<»-<";_\- 3 ) 

n=l fc=0 V ' 



[2(n - 1) + k- 1 + 2r - 2(n - l)]\ 2 r- n+ i( \M jjf ||^-«+2 



(r - n)![2(n - 1) + k - 1]! KN - 1 

< \xrr T TwjjW - (2n + fc ~ 3)! (tr + k-iy. 

~ n^tt^o k\(n-k)\ (n + fc-1)! (r-n)![2(n-l) + Jfe-l]I ' V ' 

We have ^, < < 3* and (r , ?t) ^^ fc „ 1]! < 3^*-i( r - „ + 2)! < 3*(r + 1)!. 

Moreover J2n=i J2k=o 1 < ( r + I) 2 thus for || JJ^\\ < 1 we get a bound on the contribution of the first terms 
to the reminder 

(23 3 3 3 ) r |A| r (r + 1) 3 H < (23 3 3 3 e 2 )V!|Ar . (120) 

Note that although the bound we have established might not appear tight, in fact it is: the rl growth of the 
reminder is not an artifact, but it is generated by the proliferation of the Wick contractions in a graph with 
loop lines. 

□ 

We have thus far given a non perturbative definition for the cumulants, eq. (|103[) . and showed that it 
is the Borel sum of its perturbative expansion in Feynman graphs. We now prove that the measure thus 
defined is properly uniformly bounded. 



Theorem 8. The perturbed Gaussian measure in eq. (|84[) is trace invariant and properly uniformly bounded. 

Proof: We use the invariance under unitary transformations. We add a fictitious integral over the 
unitary group U (N) , i.e. we write 

W(J, J; X,N) = [ [dU] W(J, J; A, N) , (121) 

JU(N) 
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which of course holds as nothing depends of U on the right hand side. Now, for all fixed U, we perform the 
change of variables of Jacobian 1, — > Wa^'U . in eq. (|103|) . The Gaussian measure is invariant under 
this change of variables, hence 



du l 



fc=0 i 1 ,i 2 ...,i fc = i r O i 



(i,i)er„ 



n>l ' '— n ' ' -■- ' — ^ ■'O 

Rearranging the terms in the product over d we have 

W(J,J-,X,N)=N°-^^-Xr-i±li ± eT( n **« 



n 

jeH[^-i(i)^(i)] 



j' 5 d-i (1) i S d (1) 



J d[i w i ]{Tn ^ mi {(T {1 \...(J (n) ) 



fe=0 i 1 ,i 2 ...,* fc =i 7 -0 



(i,j)er« 



fc=0 i 1 ,« 2 ...,i fe =i r 

fe 



(i,j)er„ 



jes[d->-e(<Q] 



. d=l 



3dl^{d) 



d=l 



(122) 



(123) 



(124) 



d=i jes[d-^(d)] 
The integral over the unitary group of a product of matrix elements is (see |23j ) 

/ [dU] J] U ldad Ul Jd ^WgfJV^r- 1 ) TJ Si^ w S adPrW , 

•^W d=l a,r d=l 

where a and t run over the permutations of k elements and Wg(AT, cr) is Weingarten's function [23], thus 

n>l ' '— n ' -- -- --. — fi J ,„-^^t- 



The external sources group into a product of traces. Following the indices we see that (3d —> &d —> P T (d) ■ ■ ■ 
thus each trace of a product of insertions reproduces a cycle in the permutation r. Denoting these cycles 
C r (r), and their length |C r (r)| (hence r writes t = C\{t) . . . C c ( t )(t)) we get 



c(r) 



n[^a.^ w =n tr [(^) | ^ T)| ] 



(125) 



d=l 



Similarly, the indices j, I follow the cycles of the permutation at; as jd l^M) ~* 3tr£{d) •••■> thus the generating 
function of the cumulants is 



W(J, J; A, N) = N 



D— 



n>l fc=0 i 1 ,i 2 ...,i k =l T JU fwier„ 



(y)er„ 



- c(r) 

/ ^(r„ )U) i®i(^ (1) ,--^ (n) ) E w eW £rr " 1 )II tr K Jjt ) |CV(T)l ] 
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n tr[ n ( n • ( i26 ) 

h=l d=l jeH[(^)^i(9)^(^) d - 1 (9)] 

with q any element in the cycle C/j(ctt). The cumulants are defined according to eq. (|86[) as the partial 
derivatives of W{J, J; A, N). It follows that the distribution is trace invariant, as the non trivial cumulants 
at order 2k write as sums over graphs B, whose connected components are the cycles over the external 
insertions JJ^. To each graph one has several possible t associated permutations. The number of cycles of 
r is the number of connected components of the graph $, p(B) = c(r). We now show that the cumulants 
are properly uniformly bounded. The Wcingarten function respects 

lim N^W S (N,a) = jf(-l)I^WI ( ^^j) , (127) 

N->oo |C S (CT)| + 1 V \Cs{(T)\ J 

bounding the traces of products of resolvent by N, and taking into account that the Gaussian integrals are 
normalized, and the integrals over the parameters u are bounded by 1, we get a bound for each graph B 
contributing to a cumulant of order 2k (using K(B) as a dustbin notation for a constant independent of TV, 
but depending on B) 

\Si(B,fi N )\ KK^N-^-VN - 1 ^^ 11 ^ E E E N- 2k+c ^ T ^ +c ^ (128) 

n=k i 1 ,i 2 ...,i fc = i 7-O cr,reS(fc) 

By lemma [SI c(£) + c(er£) + c(r) + c(<jt -1 ) < 2 + 2k and taking into account that £ is a cyclic permutation, 
c(£) = 1, and that the sums over r and a do not depend on n, we get a bound 

K{B)ND - mD -x)- m g |A|" ( ^*" 3) ' , , (129) 
^ (n-/c)!(n + A;- 1)! 

for some constant K(B), and the sum over n converges by the usual bounds. 

□ 



B Other scalings of the cumulants 

A natural question is to what extent the results presented in this paper can be generalized for different 
scalings scalings of the cumulants. As already mentioned the scaling N ^ 1 of the Gaussian is the unique 
scaling which leads to convergent expectations for all invariants, not only for subclasses of invariants. 

An interesting question is what happens if one allows the scaling of the cumulants to depend on finer 
details of the associated graphs. Of course if this extra scaling suppresses some of the cumulants the results 
hold. The interesting question is how much these scaling can be boosted, while still having a large N limit 
(universal or not). One particular scaling one can consider is to boost each invariant by a factor iV°' 8 ^ 



^2k 



p(B) 

-J Hl , % ■■■%,}= E N- 2(D - 1 ^+ D -^+ n ^ K(B, N) J] 6^ , (130) 

k(B)=k 

with f2(£>) its convergence order (note that the convergence order, like the degree, factors over the connected 
components of the graph fl(B) = J2 P ^(&p))- The expectation of an observable writes again as a sum over 
doubled graphs Q, 

iV- 1+0 ( B V(Tr e (T,T)) = Y, \{ K ^)) 

jy- (p-i)! "(g)+ (p-D! min S'\£0 = B"(S')+£„,pmin Sp ( )\£0 = 8p(c0 w(S P (<*))-£> J2 a (p(s(a)-l) ) (131) 
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where we have expressed the convergence orders £l(B) and fl(B p (a)) as 

= ^ g ^^ uiSpi a)) - J^^BM)) , d32) 

with Q' and Q p (a) covering graphs of B and B p (a). Again the contribution of Q is dominant only if all the 
cumulants have a unique connected component p(B{a)) = 1. that is B(a) = Bi(a), which will be the case 
we consider from now on. Let us denote the total scaling with N in eq. (|131[) 

A(£) = - 2 ttj u{G)+ ._. 2 ... min w(Q') + V mm w(S(a)) 

= ^F w (C?) + C|a|- sup sup ^F 0i (C?(«)). (133) 

If i3(a) are all dipoles B^ 2 \ corresponding to Q = g min U £0 rg m i n \ £>( 2 ' for the minimal covering graphs 
gmm Q £ ^ we Q^tajn A(^ min U£0(gmin)) = 0, as all Q(a) are of degree 0. This reproduces the usual Gaussian 
evaluation. For the other doubled graphs Q there are three possible scenarios 

• for all g D B, G + G min U £ o (gmi „) 6( 2 ), A(G) < O. In this case the model admits a Gaussian large N 
limit. 

• for all g D B, g ^ g nlin U £0{gmin) B^\ A(Q) < 0, and there exists g ^ g min U £ o (gm ,„ ) e( 2 ) with A(g) = 0. 
In this case the model admits a large N limit which is not Gaussian. 

• there exists g Z) B with A(g) > 0. In this case the model does not admit a large N limit. 

The example of figure [15] shows that we are not in the first case. It consists in an observable and a 

„ 2 



Figure 15: A doubled graph scaling like the Gaussian contribution. 



cumulant in D = 3. The minimal graphs for both the observable and the cumulant have 6 faces of colors 
Oi (hence degree w(C/(a) mm ) = u(g mm ) = 3). The doubled graph has 9 faces of colors Oi (that is degree 
cj(Q) = 6) thus A(g) = 0. 

It is for now an open question to discern in which of the remaining two cases we are. One can show 
using a Cauchy-Schwarz inequality that for any observable with 2k(B) vertices < Dk(B) for all 

doubled graphs g D B, and that the bound is saturated (by the case in which one has only one cumulant 
whose associated graph is the mirror image of the observable B, as it is the case in figure ITS]) . It follows that 
the model admits a large N limit if and only if, for all connected graphs B, one has 

m „J2 F « (S )>^Sl±^. (134) 
g\£°=B A r l 2 

i 

While we have not been able to find any example in which this inequality is violated, we have not been 
able to prove it either. 
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Should this inequality hold, one would get a non Gaussian large N limit, which of course would be very 
interesting. Note however that if the models admits a large N limit with these scalings, then the leading 
order is rather non trivial. The example in figure [15] shows that at leading order one gets contributions 
from graphs which do not correspond to manifolds. Also, it is not clear (and it does seem unlikely) that the 
leading order graphs form a summable family like the planar or the melonic graphs. 
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